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Abstract 

This paper concerns the Cauchy problem of the two-dimensional (2D) nonhomogeneous incompressible 
nematic liquid crystal flows on the whole space R 2 with vacuum as far field density. It is proved that the 
2D nonhomogeneous incompressible nematic liquid crystal flows admits a unique global strong solution 
provided the initial data density and the gradient of orientation decay not too slow at infinity, and the 
initial orientation satisfies a geometric condition (see 0). In particular, the initial data can be arbitrarily 
large and the initial density may contain vacuum states and even have compact support. Furthermore, 
the large time behavior of the solution is also obtained. 
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1 Introduction 


Liquid crystals can form and remain in an intermediate phase of matter between liquids and solids. When 
a solid melts, if the energy gain is enough to overcome the positional order but the shape of the molecules 
prevents the immediate collapse of orientational order, liquid crystals are formed. The lack of positional order 
is a shared property of liquid crystals and liquids; on the other hand, liquid crystals are anisotropic (like 
solids). The nematic liquid crystals exhibit longrange ordering in the sense that their rigid rod-like molecules 
arrange themselves with their long axes parallel to each other. Their molecules float around as in a liquid, 
but have the tendency to align along a preferred direction due to their orientation. The hydrodynamic theory 
of the nematic liquid crystals was first derived by Ericksen and Leslie during the period of 1958 through 1968 
(see mm )- A brief account of the Ericksen-Leslie theory on nematic liquid crystal flows and the derivations 
of several approximate systems can be found in the appendix of m- For more details on the hydrodynamic 
continuum theory of liquid crystals, we refer the readers to the book of Stewart 1351 . 
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In this paper, we investigate the global existence of solutions to the following two dimensional (2D) 
simplified version of nematic liquid crystal flow in the whole space R 2 , which describes the motion of a 
nonhomogeneous incompressible flow of nematic liquid crystals: (see, e.g., 125 1 126 ] !: 

g t + div(gu) = 0 , 

gut + gu ■ Vu — is A u + VP = —A div(Vd © Vd), 

( 1 - 1 ) 

dt + (u • V)d = 7(A d + |Vd| 2 d), 
divw = 0, |d| = 1, 

where g{x,t ) : R 2 x (0,+oo) —> R is the density, u(x,t) : R 2 x (0,+oo) —>■ R 2 is the unknown velocity field 
of the flow, P{x,t ) : R 2 x (0,+oo) —>• R is the scalar pressure and d : R 2 x (0,+oo) —>• § 2 , the unit sphere 
in R 3 , is the unknown (averaged) macroscopic/continuum molecule orientation of the nematic liquid crystal 
flow, divw = 0 represents the incompressible condition, and is, A and 7 are positive numbers associated to 
the properties of the material: is is the kinematic viscosity, A is the competition between kinetic energy and 
potential energy, and 7 is the microscopic elastic relaxation time for the molecular orientation field. The 
notation Vd© Vd denotes the 2x2 matrix whose th entry is given by did ■ djd (1 < i,j < 2). 

We consider the Cauchy problem for m with the initial conditions 

g(x, 0) = go(x), u(x, 0) = uo(x), d(x, 0) = do(x), |do(a;)| = 1, in R 2 (1.2) 

for given initial data go, uq, and do- Since the concrete values of is, A and 7 do not play a special role in our 
discussion, for simplicity, we assume that they are all equal to one throughout this paper. 

The above system HU E3 is a simplified version of the Ericksen-Leslie model | 8 lU 8 l , but it still retains 
most important mathematical structures as well as most of the essential difficulties of the original Ericksen- 
Leslie model. Mathematically, system ([II]) CE3 is a strongly coupled system between the nonhomogeneous 
incompressible Navier-Stokes equations (see, e.g., [24ll32lf85]. l40l 1 and the transported heat flows of harmonic 
map (see, e.g., [2] fT6ll41] ). and thus, its mathematical analysis is full of challenges. 

There is a huge literature on the studies about global existence and behaviors of solutions to (11.11) . The 
important progress on the global existence of strong or weak solutions in multi-dimension has been made by 
many authors, refer to p0i9l|2Tl|22j|42] and references therein. However, since the system HU contains the 
incompressible Navier-Stokes equations as a subsystem, one cannot expect, in general, any better results than 
those for the Navier-Stokes equations. In the homogeneous case, i.e., g = constant, Lin-Lin-Wang in [26] 
established that there exists global Leray-Hopf type weak solutions to the initial boundary value problem for 
system in]) diu on bounded domains in two space dimensions (see also m- The uniqueness of such weak 
solutions is proved by Lin-Wang [3D], see also Xu-Zhang [33] for related works. When the space dimension 
is three, Lin-Wang m obtained the existence of global weak solutions very recently when the initial data 
(uo,do) M 2 x H 1 with the initial director field do maps to the upper hemisphere §1. The global existence 
of weak solutions to homogeneous system (O) ([H with general initial data in dimension three is still an 
open problem. There are also some studies on the homogeneous incompressible nematic liquid crystal flows, 
we refer the readers to [H[23|28l[33] and references therein. 

In the non-homogeneous case, i.e., the density dependent case, when the initial data away from vacuum, 
recently, Li [21 established the global existence of strong and weak solutions to the two-dimensional system 
( 11 . 11 ) ( 11 . 21 ) provided that the initial orientation do = (doi, do 2 , ^ 03 ) satisfies a geometric condition 

do 3 > £0 for some positive £0 > 0. (1.3) 

In the presence of vacuum, if the initial data is small (in some sense) and satisfies the following compatibility 
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conditions 


~liAu 0 + VP 0 + A div(Vd 0 © Vd 0 ) = £>o 9o (1-4) 

in a bounded smooth domain Q C M 2 , and {Po,go) G x L 2 (S1), Wen-Ding [32] obtained the global 

existence and uniqueness of the strong solutions to the system (PD PD . see also mm for related works. 
It should be emphasized the possible appearance of vacuum is one of the main difficulties, which indeed leads 
to the singular behaviors of solutions in the presence of vacuum, such as the finite time blow-up of smooth 
solutions un¬ 
it is not known in general about the existence of global strong solutions to the problem PD PD in 
two-dimension without the geometric condition m or the compatibility condition imposed on the 

initial data. This problem is rather interesting and hard to investigate. Indeed, it should be noted that the 
previous studies on the heat flow of a harmonic map [ 2 ] indicate that the strong solution of a harmonic map 
can be blow-up in finite time. In our case, since the system contains the heat flow of a harmonic map as 
a subsystem, we cannot expect that m have a global strong solution with general initial data. This makes 
the analysis rather delicate and difficult. 

It should be noticed that when d is a constant vector and Ml = 1, the system (11.11) reduces to the 
nonhomogeneous incompressible Navier-Stokes equations, which have been discussed in numerous studies 
[H151 1 15115^1 136) and so on. In the case when go is bounded away from zero, Kazhikov m established the 
global existence of weak solutions. Danchin jJJ] obtained the existence of local and unique strong solutions in 
the so-called critical Besov space. When the initial data may contain vacuum states, Lions [32] obtained the 
the global existence of weak solutions. Choe-Kim [3] proposed a compatibility condition and established the 
local existence of strong solutions under some suitable smallness conditions. The global existence of strong 
solutions on bounded domains were established by Huang-Wang M- Recently, inspired by Li-Liang ns. 
Lii-Shi-Zhong [36] established the global existence of strong solutions to the 2D Cauchy problem of the 
incompressible Navier-Stokes equations on the whole space R 2 with vacuum as far field density. However, the 
global existence of strong solution with large data to the 2D Cauchy problem m dm with vacuum as far 
field density is still open. In fact, this is the main aim of the present paper. 

Before formulating our main result, let us first define precisely what we mean by strong solutions. 

Definition 1.1 If all derivatives involved in (11.11) for (p, u, P, d) are regular distributions, and equations (11.11) 
hold almost everywhere in R 2 x (0,T), then (g,u,P,d) is called a strong solution to (11.11) . 


Without loss of generality, we assume that the initial density go satisfies 



= 1 , 


which implies that there exists a positive constant Nq such that 



£>odx > - 



1 

2 ’ 


where Br — {x G R 2 ||x| < R} for all R > 0. 
Now, we state our main result as follows: 


(1.5) 


( 1 . 6 ) 


Theorem 1.2 For constants q > 2, a > 1, assume that the initial data (go,uo,do) satisfies (11.31) . (11.51) . (11.61) . 
and 


go > 0, g 0 x a sL 1 (R 2 )ClL 1 (R 2 )nW 1 'P(R 2 ), p 0 ! a 0 ei 2 (R 2 ), V « 0 £ i 2 (R 2 ), 


divuo = 0, do £ L 2 (R 2 ), Vd 0 x? G L 


V 2 d 0 G L 2 


Mo| — i, 


(1.7) 
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where 


x = (e + |x | 2 ) 2 ln 2 (e + |x| 2 ). 


T/ien t/ie Cauchy problem Oil (fL 2 l) has a unique global strong solution ( g,u,P,d ) satisfying that for any 
0 < T < oo , 

0 < 0 G C([0, T]; L\R 2 ) n R^R 2 ) n W^QR 2 )), 

gx a G T°°(0, T; T^R 2 ) n R^R 2 ) n W^R 2 )) 

g^u, Vu, ti S7u, ti giu,t^S/P,tS7P, t^S7 2 u, tV 2 u G T°°(0, T ; T 2 (R 2 )), 

S/d, Vdif , V 2 d,tiv 2 d,tiVdt,tiv 3 d G L°°(0, T; T 2 (R 2 )), 

< VMGL 2 (0,T; J ff 1 (R 2 ))nL^(0,r ; W^ 1 ’«(M 2 )), (1.8) 

VP G L 2 (0, T; T 2 (R 2 )) n L^r (0, T; T«(R 2 )), 

V 2 d G T 2 (0,T;P 1 (R 2 )), Vdt, V 2 ctet g P 2 (R 2 x (0,T)), 
g?u u t?S7u t ,t?S7d t ,t?S7 2 d t G L 2 (R 2 x (0,T)), 
t*Vu G T 2 (0,T; W' 1 >«(R 2 )), 

< 

and 

inf [ g(x,t)dx > -7 (1.9) 

0 <*<tJ Bni ~ 4 

for some positive constant N± depending only on ||£>g uo||l 2 (k 2 ); -Wb and T. Moreover, the solution ( g,u,P,d ) 
has the following temporal decay rates, i.e., for all t > 1 , 


| || <? 2 (', 0 II L 2 (R 2 ) + IMt(’ji)||fl' 1 (K 2 ) + ||Vu(-, t)||i2( K 2) + ||V 2 d(-, f)|| L 2( R 2) < Ct 2, 
\||V 2 u(-,f)|| L 2 (R2) + |||Vd||V 2 d|(.,<)|| L 2 (R 2 ) + ||VP(-,t)|| i 2 (R 2) < Ct-\ 

where C depends only on initial datums. 


Remark 1.3 When d is a constant vector and |d| = 1, the system (11.11) turns to be the nonhomogeneous 
incompressible Navier-Stokes equations, and Theorem 11.0 is similar to the results of 136V . Roughly speaking, 
we generalize the results of fSffl to the incompressible nematic liquid crystal flows. 

Remark 1.4 Our Theorem ! 1. 21 holds for arbitrarily large data which is in sharp contrast to ffM21l\f2} where the 
smallness conditions on the initial data is needed in order to obtain the global existence of strong solutions. 
Moreover, the compatibility condition CEO) on the initial data is also needed in flA \21\\lfl^ . It seems more 
involved to show the global existence of strong solutions with general initial data. This is the main reason for 
us to add an additional geometric condition m- 

Remark 1.5 Compared with 120/ . there is no need to impose the absence of vacuum for the initial density. 
Furthermore, it should be pointed out that the large time asymptotic decay with rates of the global strong 
solution in (11.101) is completely new for the 2D nonhomogeneous nematic liquid crystal flows. 

We now make some comments on the analysis of the present paper. Using some key ideas due to |241I35I . 
where the authors deal with the 2D incompressible Navier-Stokes and MHD equations, respectively, we first 
establish that if (go, uq, do) satisfies (11.51) (11.71) . then there exists a small T 0 > 0 such that the Cauchy problem 
Ob ED admits a unique strong solution ( g,u,P,d ) in R 2 x (0, To] satisfying (11.81) and (11.91) (see Theorem 
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ED- Thus, to prove Theorem II.21 we only need to give some global a priori estimates on the strong solutions 
to system (friD lira in suitable higher norms. It should be pointed out that the crucial techniques of proofs 
in EEDl cannot be adapted to the situation treated here. One reason is that when fl C R 2 becomes unbounded, 
the standard Sobolev embedding inequality is critical, and it seems difficult to bound the L p -norm (p > 2) 
of the velocity u just in terms of ||{>*m||l 2 ( H 2 \ and || Vzi||i 2 ( R 2 ). Moreover, compared with [35][36], for system 
(ED ED treated here, the strong coupling terms and strong nonlinear terms, such as div(Vd © Vd), u • Vd 
and |Vd| 2 d, will bring out some new difficulties. 

To overcome these difficulties mentioned above, some new ideas are needed. To deal with the difficulty 
caused by the lack of the Sobolev inequality, we observe that, in equations ED, the velocity u is always 
accompanied by g. Motivated by mmm, by introducing a weighted function to the density, as well as the 
Hardy-type inequality in |32] by Lions, the (r > 2, a > 0 and H = Br or R 2 ) is controlled in terms 

of || e ^ u \\l 2 {Q) and ||(see Lemma [2~4ll . After some spatial estimates on Vd (i.e., Vdx?, see (13.511 ). 
and suitable a priori estimates, we then construct approximate solutions to ED, that is, for density strictly 
away from vacuum initially, we consider a initial boundary value problem of ED in any bounded ball Br 
with radius R > 0. Combining all these ideas stated above with those due to [111135], we derive some desired 
bounds on the gradients of the velocity and the spatial weighted ones on both the density and its gradients 
where all these bounds are independent of both the radius of the balls Br and the lower bound of the initial 
density, and then obtain the local existence and uniqueness of solution (see Subsection 3.2). 

Base on the local existence result (see Theorem l3.1D , we try to give some a priori estimates which are needed 
to obtain the global existence of strong solutions. When we derive the estimates on the L°°(0, T; L 2 (R 2 ))- 
norm of || Vu||£ 2 ( R 2 ) and || V 2 d||i 2 ( R 2 ). On the one hand, motivated by |2C)][36], we use material derivatives 
u = ut+u-Vu instead of the usual ut , and use some facts on Hardy and BMO spaces (see Lemma l2?6ll to bound 
the key term f R2 |P||Vu| 2 da; (see the estimates of h of (14.51) 1. On the other hand, the usual L 2 (R 2 x (0,T))- 
norm of V dt cannot be directly estimated due to the strong coupled term u ■ S7d and the strong nonlinear 
term |Vd| 2 d. Motivated by [37], multiplying (II. ll) 2 by A Vd instead of the usual Vd t , and the nonlinear terms 
u ■ X7d and \Vd\ 2 d can be controlled after integration by parts (see (14.91) 1. Using the structure of the 2D 
heat flows of harmonic maps, we multiply (13.81) by VdA| Vd| 2 and thus obtain some useful a priori estimates 
on |||Vd||V 2 d||| i 2 ( R 2 ) and |||Vc?||AVd|||z, 2 ( R 2 ) (see (14.151) 1. which are crucial in deriving the time-independent 
estimates on both the L°°(0, T; L 2 (R 2 ))-norm of tig^u and the L 2 (R 2 x (0,T))-norm of t^Vii (see (14.111) 1. 
Next, after some careful analysis, we derive the desired L 1 (0, T ; L°°(R 2 )) bound of the gradient of the velocity 
Vu (see (14.281) 1. which in particular implies the bound on the L°°(0, T; L 9 (R 2 ))-norm (q > 2) of the gradient 
of the density. Moreover, some useful spatial weighted estimates on g, Vd, V 2 d are derived (see Lemma [ED- 
With the a priori estimates stated above in hand, we can estimate the higher order derivatives of the solution 
(p, u, P, d) (see (14.331) 1 by using the similar arguments as those in [3B]E3 to obtain the desired results. 

The remaining parts of the present paper are organized as follows. In Section 2, we shall give some 
elementary facts and inequalities which will be needed in later analysis. In Section 3, we prove the local 
existence and uniqueness of solutions to the Cauchy problems of ED ED- Finally, in Section 4, we first 
establish some a priori estimates on strong solutions, then the proof of Theorem [L2] is given. Throughout the 
paper, we denote by C the positive constant, which may depend on a, £o and the initial data, and its value 
may change from line to line. We also C(a , /3, ■ • •) to emphasize that the constant C depends on a, /3, ■ ■ ■. 

2 Preliminaries 

In this section, we shall give some known results and elementary inequalities which will be used frequently 
later. We first list the following local existence theory on bounded balls, where the initial density is strictly 
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away from vacuum, whose proof can be shown by similar arguments as in [1311211142] . 
Lemma 2.1 For R > 0 and B R = {i£ K. 2 1 |ct| < R}, assume that (po,Uo,do) satisfies 


go eH\B R )nL°°{B R ), u 0 &Ho(B R )nH\B R ), d 0 £H 2 (B R ), inf p o (x)>0, div Uo = 0, |d 0 | = 1. (2.1) 

xGB R 

Then there exists a small time T R > 0 such that the equations ED with the following initial-boundary value 
conditions 


( 2 . 2 ) 


(g,u,d)(x,0) = (g 0 ,u 0 ,d 0 ), x £ B R , 
u(x, t) = 0 ^ d(x, t) = 0, x £ dB R , t > 0, 

has a unique solution (g,u,P,d) on B R x (0 ,T R \ satisfying 

' g £ C((0, T r ]-H' 2 (B r )) n L°°(B r x (0, Tfl]), g t G L°°(0, T; L 2 (B R )), 
u G L°°(0,T R -,Hi(B R )nH 2 (B R )), u t G L 2 (0,T R ; H^B R )), 

P£L°°{0,T r -H\B r )), 

d G L°°(0,T R ; H 3 (B R )), d t £L°°(0,T R -H 1 (B R )). 

Here, v = (with \x\ = R) is the outer normal vector on dB R . 

Next, for either 12 = R 2 or 12 = B R with R > 1 , the following weighted L p -bounds for elements of the 
Hilbert space D l ^ 2 (Ll) = {v £ Hl oc (Ll)\\7u £ L 2 (12)} can be found in Theorem B.l of [32]. 


(2.3) 


Lemma 2.2 For m £ [2, oo) and 9 £ ( 1+ 2 m , oo), there exists a positive constant C such that for either 12 = R 2 
or 12 = B r with R> 1, and for any v £ D 1,2 ( 12), 

(X (ln(e + ^ 2)yBdx ) m - + C7||v«iu a(n) . 

A useful consequence of Lemma [2~2l is the following weighted bounds for elements of Z2 1 , 2 (12), which have 
been proved in [191124] . It will play a crucial role in our following analysis. 


Lemma 2.3 Let x be as in Theorem \1.2\ and 12 as in Lemma I£.21 Assume that g £ L 1 ( 12) n L°°(12) is a 
non-negative function such that 



gdx > Mi, 


||e||i, 1 (n)ni°°(n) < M2, 


for positive constants Mi, M 2 , and Ni > 1 with B jv x C 12. Then for e > 0 and r/ > 0, there is a positive 
constant C depending only on e,rj, Mi, M 2 , and Ni such that every v £ Z2 1 , 2 (12) satisfies 


\\vx ’ 7 ll i ^±£ (n) < C\\e 2 v\\ L 2 (Q) + C||Vv|| i 2 (n ) 


with rj = min{l, 77 }. 


We shall still need to the following lemma obtained by Liang [2Tj . 

Lemma 2.4 Let the assumptions in Lemm,a \2.S\ holds. Suppose in addition that gx a £ L 1 ( 12) with a > 1. 
Then there is a constant C depending on Mi, M 2 , Ni,a, r, ||p||i,°°(n) o,nd ||pal a ||Li(n) such that for a > r ( 2 + a ) 
and r > 2 

< C(ll£’^'y||L 2 (f2) + ||Vu|| L 2 (a) ). (2.4) 
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Let us recall the following L p -estimate for elliptic systems, whose proof is similar to that of [3] (see Lemma 
12). It is also a direct result of the combination of the well-known elliptic theory [L and a standard scaling 
procedure. 


Lemma 2.5 For p > 1 and k > 0, there exists a positive constant Co depending only on p and k such that 


l|V fc+2 v|| lp(b r ) < C'o||Au|| V yfe, P ( BR ) 


for every v £ W k+2,P (B R ) satisfying either 


v ■ x = 0, rotu = 0, on dB R 


(2.5) 


or 


n • Vu = 0, on 8 Br. 


Let 'H 1 (K 2 ) and BMO( K 2 ) stand for the usual Hardy and BMO spaces (see [S]), we end this section by 
the following useful lemma, whose proofs can be found in [5] . 

Lemma 2.6 (i) There is a positive constant C such that 


\\E ■ -®||'H 1 (h 2 ) < C||I1'IIl 2 (r 2 )I|I3||l 2 (r 2 ) 


for all E £ L 2 (R 2 ) and B £ L 2 (R 2 ) with 


div E = 0, V ± - B = 0 in V(M. 2 ). 


(ii) There is a positive constant C such that for all v £ D 1 , 2 (K 2 ), it holds 


v\\bmo(wl 2 ) < C||Vz;||z ; 2( R 2). 


3 Local existence and uniqueness of solutions 

In this section, we shall prove the following local existence and uniqueness of strong solutions to the Cauchy 
problem of (ED ED - 

Theorem 3.1 Assume that (ffoiUo, do) satisfies (11.51) (11.71) . Then there exist a small positive time To > 0 
and a unique strong solution ( g,u,P,d ) to the Cauchy problem of system (11.11) (11.21) in R 2 x (0,To] satisfying 
(11.81) and (11.91) with T = Tq. 

Remark 3.2 We remark that in Theorem no we do not need the condition ED- From mm, it is easy 
to see that if we replace x in Theorem \1.2\ by 


(e + |cc| 2 ) 2 In 1-1 " 170 (e + |cc| 2 ), 


with rjo > 0, then the result of Theorem \3.1\ still holds. 


3.1 A priori estimates 

Throughout this subsection, for p £ [1, oo] and k > 0, we denote 



i ’ p = W 1 ' p {B r ), H k = W k ' 2 
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for simplicity. Denote by || • ||x the norm of the X(i3/{)-functions. Moreover, for R > 4Nq > 4, in addition 
( 12 . 211 . assume that (go,uo,do) satisfies 

\ < f £>o(a;)dx < f p 0 (^)da; < 1. (3.1) 

^ J Bn 0 J Bb 

Then, from Lemma |2.II we find that there exist some Tr > 0 such that the initial-boundary value problem 
(ini) and (12.21) has a unique strong solution ( g,u,P,d ) on Br x [0,Tr] satisfying (12.31) . 

For x, a and q as in Theorem fi~2l the main goal of this subsection is to derive the following key a priori 
estimates on <h(f) defined by 

$(f) = 1 + \\giu\\- L 2 + ||Vu||£ a + ||Vdi 2 1|| 2 + ||V“d||| 2 + ||eS a ||L 1 nff 1 nw 1 -'>) 

which are needed for the proof of Theorem 13.11 


Proposition 3.3 Assume that (go,uo,do) satisfies ( 12 . 11 ) and (13.11) . Let ( g,u,P,d ) be the solution to the 
initial-boundary value problem CCD and ( 12 ^ 1 ) on Br x [0 ,Tr] obtained by Lemma \2.1\ Then there exist 
positive constants Tq and M both depending only on a,q : No, and Eq such that 


sup (&(t) + t,\\g 2 ut\\ 2 L 2 + t\\dt \\ 2 h i + f || V 2 u ||22 + t\\ VP|j £ 2 + t|| V 3 d||£ 2 'j 

Ct<T„ \ / 

(ll^foll ! 2 + liv 2 u||l 2 + livdtii^x + ||v 3 d ||| 2 + ||v 2 dTt||| 2 ) dt 


0 <t<T 0 
rT 0 


CTO / g+l l+g 

||V 2 u||^ + ||VP|| L l +t\\V 2 u\\l q +t\\\7P\\l q ) dt 


/o 


[ ° t (l|Vut || 2 2 + || Vdtllj? 1 ) dt < M, 

Jo 


(3.2) 


where 


Eo — 1 + ||£>o«o||z , 2 + ||Vmo||l 2 + HVdo ^ 2 ||l 2 + ||V 2 do||i 2 + || 0 oa ;a ||L 1 nBB 1 nw 1 ’P- 


To prove Proposition 13.31 whose proof will be postponed to the end of this subsection, we need to give the 
following useful estimates. 

Lemma 3.4 Under the assumptions of Provosition 1 3.3\. let ( g,u,P,d ) be a smooth solution to the initial 
boundary value problem (EH) and (12.2p . Then we have for all t >0, 

sup (\\q 2 u \\\ 2 ^-\\^d\\\^ ^ f \\\7u\\ 2 L 2 + \\Ad + \\7d\ 2 d\\ 2 L2 dT < \\q§uo \\ 2 L 2 + \\\7do\\ 2 L 2, (3.3) 

0 <r<t ' ' ./n 


sup 

0 <T<t 


L^rL°° + 


HVdlliO + J* ||V 2 d || 2 2 dr < Cexp \c $(r)dr| , (3.4) 


and 


sup (||Vdit||2 2 ) + f ||V 2 dxt||| 2 rf T < Cexpj f $(r)drl . 
0 <T<t J 0 Uo J 


(3.5) 


Proof. (13.31) is the standard energy inequality, see, e.g., |16ll26] . To prove (13.41) . we first notice that from the 
continuity equation CEDi and the divergence free condition (EH> 4 , it is easy to see (cf, [32]) that 


\\e(t)\\LP = IMIlp, for all p G [l,oo] and t> 0. 

Next, multiplying (EH 3 by —Ad, and integrating over Br , it follows that 

i^H Vd|| i a +| 11 v2 ^Hi a < / |Vu||Vd||Vd|dx + J \X7d\ 2 \V 2 d\dx + C\\Vd \\ 2 L 2 


(3.6) 


(3.7) 


<||Vu|| i2 ||Vd|| 2 4 + ||Vd|| 2 4 ||V 2 d|| i2 + C||Vd|| 2 2 
<\\\^ 2 d\\h + C{ 1 + ||V U ||| 2 + ||Vd|| 2 2 )||Vd|| 2 2 
<\\\V 2 d\\l 2 +C<S>mVd\\l 2 , 

where we have used the condition |d| = 1, and the following inequalities 


n 2 « 2 « r\ 

/ |Ad| 2 dx = ^ / d 2 d ■ djddx = — ^ / did ■ did 2 ddx (notice that — d = 0 on BBr.) 


i,j =i 


i,j —1 


= — ^ / <%d • djdidvjdS + ^ / |<9i<9,d| 2 dx (^ is defined in Lemma l2Tl f 


*.j=i 




hi=i' 


2 2 
=||V 2 d||| 2 + V [ did ■ djddiVjdS = ||V 2 d||| 2 + 4 V f dd-djdSijdS 

i,j =i ^ i,i=i 


>11 V 2 <i||22 — IIVd||| 2 ( aBR ) > ||V 2 d||| 2 — C||Vd||^i^ ^ (see, e.g., |T7] for trace embedding) 

1 , 


— I|V a||i 2 - C||Vd|| i 2 ||V 2 d || i2 > -||V 2 d ||£ 2 - C7||Vd||£ 2 . 

Here S,j = 1 if i = j (or = 0 if i ^ j). The estimate (13.71) together with Gronwall’s inequality ensures 

sup (||Vd|| 2 2 ) + / ||V 2 d|| 2 2 dr < exp ( C f <f>(r)drj . 

0 <T<t Jo L Jo J 

This combined with (13.61) leads to (13.41) . 

In what follows, we are in a position to prove (13.51) . By applying V on (11.11) .,. we have 

Vd t - AVd = — V(m • Vd) + V(|Vd| 2 d). 

Multiplying (13.81) with S7dx a and integrating by parts over Br yield 

~||Vdsc*|| 2 a + ||V 2 dx^||| 2 <cJ\\7d\\V 2 d\Vx a dx+J\Vu\\Vd\ 2 x a dx+J\u\\\7d\ 2 Vx a dx 

+ J |Vd| 2 |V 2 d|i a dx+ J\Vd\ 3 Vx a dx+C\\Vdx%\\ 2 L2 
=h + I 2 + h + h + h +C\\ Vdx% ||| 2 , 
where we have used the following inequalities 

— J VAd ■ \7dx a dx = — ^ J did 2 d ■ didx a dx 

i,j= i J 

2 2 2 
= — ^ j didjd ■ didx a VjdS + ^ J \didjd\ 2 x a dx + ^ J didjd ■ diddjX a dx 

■i,j—l i,j =1 i,j=1 

2 da;^||| 2 + f didjd ■ diddjX a dx + T'' f did ■ didx a diVjdS 


(3.8) 


(3.9) 


=IIV 


*,i=i 

2 


*d=i 


= ||V 2 da;21|| 2 + f didjd-diddjX a dx +-=; / did ■ didx a 5ijdS 

i J >' i ~ 1 JdB R 

>||V 2 di^ ||| 2 — C J\V 2 d\\Vd\Vx a dx~\\Vdx a \\ 2 L2{dBR) 
>\\V 2 dx%\\ 2 L2 -C J \V 2 d\\Vd\Vx a dx - C\\Vdx a \\ 2 H i 


?{Br) 
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>i||V 2 ctef ||| 2 -C J |V 2 d||Vd|V5 a d:r - C'||Vcfe a || 2 2 . 

We can bound each term on the right-hand side of (13.91) as follows 

h<J\Vd\\V 2 d\x a dx < ^||V 2 dx®|||2 + C'||Vd®*||£ 2 ; 

h <||V«|| ia ||Vd®T||i 4 < \\\7u\\ L 2\\Vdx%\\ L 2{\\Vdx%\\ L 2 + ||V(Vdi^)|| L 2) 
<\\Vu\\ L 2\\Vdx%\\ L 2(\\Vdx%\\ L 2 + ||V 2 dxt|| i2 + ||VdVsS|| L 0 

<^l|V 2 d.Tt||| 2 + C(l+||Vu||i 2 )||V^f||| 2 <^||V 2 ^f|| 2 2 +Cd>(t)||Vdif|| 2 2 ; 

I3 < J \u\\S 7 d\x a ~i&x <\\S7dx%\\ L i\\S7dx%\\ L 2\\ux~i\\ L i 

<C\\Vdx%\\ 2 Li + C(||e*u||£ a + ||Vu||£ 2 )||Vdir*||£ a < l||V 2 da;t|| 2 2 + C$(t)||Vda;t|| 2 2 ; 
h <||V 2 d^|| i2 ||Vdst|| i4 ||Vd|| i 4 < l||V 2 dic*||£ 2 + ||Vd|| 2 4 ||Vdif ||| 4 

<^l|V 2 dif ||| 2 + (||Vd|| 2 2 + ||V 2 d|| 2 2 )||Vdit|| i2 (||V 2 dxt|| i2 + ||VdVic*M 
<^\\V 2 dx^\l,+Cm\\Vdx%\\l 2 ; 

h<j\Vd\ 3 x a dx< ||W|| L2 ||Vd^|| 2 4 < ^||V 2 ^^|| 2 2 +C$(i)||V^||| 2 , 

due to the energy inequality (13.31) . Holder’s and Gagliardro-Nirenberg inequalities. Hence, inserting the 
estimates of Ii(i = 1,2, • • • ,5) into (13.91) . we obtain (13.51) after by using Gronwall’s inequality. This completes 
the proof of Lemma 13.41 □ 


Lemma 3.5 Under the assumptions of Proposition \3.S\. let ( g,u,P,d ) be a smooth solution to the initial 
boundary value problem (EH) and (|2.2|) . Then there exists aTi = T{Nq,Eq) > 0 such that for all t G (0, Ti], 


sup (||pa; Q || L i + ||Vit ||| 2 + ||V 2 d|| 2 2 ) 

0<t<£ 

+ / (||^u T ||| 2 + ||Vd T ||| 2 + ||AW||| 2 ) dr<C eX p $(r)d? 
j ||V 2 M|| 2 2 dr < Cexp |c J $(r)dr| . 

Proof. First, for N > 1, let ifiN G C£°(Bn) satisfy 


N 


0 < (fiN < 1 , <Pn(x) = 


1, if |x| < 2 , 
0, if |ie| > N, 


|V>iv| < CN-* for k G N. 


It follows from 0! and the energy inequality (13.31) that 


-jj / Qip 2 N 0 dx = I gu- X7tp2N 0 dx > -CN 0 


-1 


gdx 


Integrating the above inequality and using (EH), it follows that 


inf 

0<*<Ti Jg' 


g\u\ 2 dx) >-C(N 0 ,E 0 ). 


1 


gdx > inf / gip 2 N 0 dx > / g 0 <p2N o dx - CTi > -, 


2 N 0 


0<t<Ti 


(3.10) 

(3.11) 

(3.12) 


(3.13) 


where T\ = min{l, -jh}. From now on, we will always assume that t < Tf in this subsection. The combination 
of (13.131) . (13.31) and Lemma IT7TTT1 ensures that for e > 0, g > 0, and every v G D 1}2 (B^ r) satisfies 


\\vx ''ll 2 +e < C{£,r[)\\g^v \\ L 2 +C(e, 77)11 Vw|| L 2 , 
L ») 


(3.14) 
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where r/ = min{ 1, y}. Using the above estimate (13.141) . after multiplying (11.11^ by x a and integrating by parts 
over Br, we have 


-jj\\Qx a \\ L i <C j g\u\x a 1 ln 2 (e + \x\ 2 )dx < C\\gx a 1+ 


4 +“ a+4 \\UX *+* r. + 4 

"lots" 


<ciip|irL“ii^ ai,:l 


L 1 \\\B 2u Wl 2 + IIVulU 2 ) 

<C(l + ||pr|| il )(l + ||Vu|| 2 i2 ) 

due to (13.31) and (13.61) . This combined with Gronwall’s inequality and (13.31) leads to 


sup || gx a Wl 1 < Cexp 

0 <T<t 


H 


(1 + l|V«||iOdr \<C. 


Now, multiplying (0) 2 by Ut and integrating by parts over Br, one has 

— ||Vu||| 2 + ||p 5 Ut || 2 2 < u J g\ u \\Vu\\u t \dx + /(VdQVdJ.V^d*. 

By using Holder’s and Gagliardo-Nirenberg inequalities, and (12.41) ensure that 

J g\u\\Vu\\u t \dx <^\\g^u t f L 2 + C J g\u\ 2 \Vu\ 2 dx < ^\\g?u t \\ 2 L 2 + C\\g?u\\ 2 L s\\\ 7 u\\\ 
^ll^lll 2 + C'||^u|||s||Vu|| 2 2 ||Vu ||^ 1 

<2 ll£ , 2 u i|li 2 + e||V 2 u || 2 2 + G(||p 2 u ||| 2 + ||Vu||| 2 )3 ||Vm || 2 2 
<-j\\e 2 u t \\ 2 L 2 + e||V 2 u || 2 2 + C'H > 3 (t), 


(3.15) 


(3.16) 


(3.17) 


where (and in what follows) we have used /3 > 1 to denote a generic constant, which may be different from 
line to line. For the second term on the right-hand side of (13.161) . we have 

J (Vd 0 Vd) • Vittda? =— J (Vd © Vd) • Vudx — J Vd* © Vd • Vuda; — J Vd 0 Vd t • Vudx 


<- 


' dt 


(Vd © Vd) • Vitdx + -||Vdt||i2 + C||Vd|| 2 4 ||Vu|| 2 4 


-Jt / (W ° vd ) • Vudx + \\\^ d tWl 2 + C'||Vd|| L2 ||Vd|| H1 ||V«|| La ||Vu|| ffl 


<- 


' dt 


(Vd © Vd) • Vudx + h\ Vdt||| a + eirV 2 u||| 2 + C'$' 3 (t). 


Inserting (13.171) and (13.181) into (13.161) . it follows that 


—i?(t) + ||p 2 u t || 2 2 < e||V 2 u|| 2 2 + -||Vd t || 2 2 + C<&P(t), 


(3.18) 


(3.19) 


where 


satisfies 


B{t)± 



(Vd © Vd) • Vrtdx 


|liv«||i 3 - cyvdHfp < B(t) < |||v^||i 2 + c'livdn^i. 

Next, from (13.8fi . it is easy to obtain that 

|(l|V 2 d|| 2 2 + i||Vd|| 2 , (aSK) ) + ||Vd t ||B ! + ||AVd||^=y'|Vdi-AVd| 2 d.T<Cy |V(u • Vd)| 2 + |V(|Vd| 2 d)| 2 dx 
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<C||Vu||| 2 ||Vd|||^ +C||| U ||V 2 d |||| 2 +C'|||Vd||V 2 d |||| 2 +C||Vd|||a 

<C(||V M ||^||Vd|b(||Vd|b + ||V 3 d|| z ,) + ||ux-^|| 2 8 ||V 2 dxt||^||V 2 d || i 4 + ||Vd|| 2 4 ||V 2 d || 2 4 + ||Vd|||e) 
<C , (||V«||i 2 ||Vd|| ia ([[Vd|| i 2 + ||V 3 d|| L 2) + \\ux~l Hia||V 2 dx 5 || i 2 (||V 2 d || i2 + ||v 3 d|| i2 ) 

+ || Vd||i 2 1| V 2 d||| 2 (|| V 2 d || i 2 + ||V 3 d|| L 2 ) + ||Vd|| 2 2 ||V 2 d||i 2 ) 

-i^o l|V3d|1 ' 2 + C||v2 ^ f Hi 2 +C* p (t), (3.20) 

where we have used Holder’s and Gagliardo-Nirenberg inequalities, (13.141) . and the following equalities 


-2 


Vd t AVddx = 


^ r a 

2 Y / didt ■ didjdujdS + -||V 2 d ||| 2 

ij=i J 9 b r at 


2 . , 

=2 Y / ■ djddiVjdS + — : ||V 2 d|| 

i,j — l < U 

2 

=|E/ ' djdSijdS + E|V 2 d|| 

U , ;_1 JOBr at 


i,j=1 JdB R 

= El|V 2 d||i2 + i||Vd||i 2(aBR) ), 


dt 


R' 


2 

L 2 


2 

L 2 


due to J^d = 0. Here, Co is the constant defined in Lemma T2. 5 1 On the other hand, by letting v = Vd, then 
it is easy to see that the condition (12.51) holds. Hence, Lemma [2.51 ensures that 


||V 3 d || L 2 <C 0 ||AVd|| L 2 . 


Inserting the above inequality into (13.201) . and then adding the resulting inequality to (13.71) . it follows that 
^(l|Vd|| 2 ffl + -^||Vd|& (aiJj0 ) +||Vd t ||^ + i ^||V 3 d || 2 2 <C||V 2 d ^||| 2 + C&(t). 

Multiplying the above inequality by C\ + 1, and then adding the resulting inequality to (13.191) . it holds that 

— (H(t)+(C'i + l)(|| VdH^i + —1| Vd|| 2 2 ( aSB ))) + ||d 2 M t || 2 2 + ||Vd t || 2 2 + ^£rl|V 3 d ||| 2 

<e|| V 2 m||| 2 + C\\V 2 dx% 11^2 + C$P(t). (3.21) 

On the other hand, notice that (g, u, P , d) satisfies the following Stokes system 

—Am + VP = — gut — gu ■ Vu + div(Vd © Vd), x £ Bn, 

div u = 0 , x £ Bn, 

u(x) = 0 , x £ dBn- 

Applying the standard L p -estimate to the above system (see, e.g., [40]) yields that for any p £ ( 1 , oo), 

||V 2 u||lp + I|VP|| LP < C\\gut\\LP + C\\gu • Vu\\lp + C*|| div(Vd © Vd)||i P , (3.22) 


from which, after using (12.41) . (13.31) . (13.61) . and Gagliardo-Nirenberg inequality, we have 

||V 2 M||i 2 + || VP || 2 2 <C ||^ t ||| 2 + C\\gu • Vu || 2 2 + C|| div(Vd © Vd )|| 2 2 

<C||e||^||^M t || 2 2 +G||^|| 2 4 ||VM || 2 4 + G||Vd|| 2 4 ||V 2 d || 2 4 
<c||^ Ut || 2 2 +( 7 (||^ u || 2 2 + ||v u ||| 2 )||v u || i 2 ||VM|| J? i+q|vd|U 2 ||v 2 d||l 2 ||v 2 d||^ 

A^=rl|V 3 d ||! 2 + -||V 2 m|| 2 2 + C\\g 2 u t \\n 2 + C{ 1 + ||Vu ||® 2 + || N 7 ^d||^ 2 ) 
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— 4C70 II V3d H^ 2 2 ll^” u lli 2 + ^'ll£ ,2 ' u t|l! 2 + C^(t). 


( 3 . 23 ) 


Substituting (13.231) into (13.211) . and choosing s suitably small, we get 

— {Bit) + (C'i+l)(||Vd||^i + —1| + 2llfi' 2u tllL 2 + ll^tlll 2 + 2 C^ H^ 3 ^Hi 2 

<C\\V 2 dx%\\ 2 L 2 +C<S> p (t). 


Integrating the above inequality with respect to time variable over (0, t), it follows from the definition of B(t), 
(13.51) and Lemma 12.51 that (13.101) holds. 

By using (13.231) and (13.101) . it is easy to see that (13.111) holds. This completes the proof of Lemma RT5l □ 


Lemma 3.6 Let (g,u,P,d) and T\ be as in Lemma Iff. 51 Then there exists a positive constant (3 > 1 such 
that for all t G (0,Ti], 

sup (r\\Q^u T \\ 2 L 2 +T\\d t \\ 2 H i S j+ J(T\\Vu T \\ 2 L 2 +T\\Vd T \\ 2 H i)dT <Cexp jcexp i^cj $ / 3 (r)dr||. (3.24) 

Proof. Differentiating (11.11) ^ with respect to time variable t gives 

gutt + gu ■ S7ut — A ut = —gt{ut + u ■ Vu) — gut ■ Vw — VPt — div(Vd © Vd)*. (3.25) 

Multiplying (13.251) by ut and integrating the resulting equality by parts over Bn, we obtain after applying 
(EDr and the divergence free condition ( 11 . 11 ) ^ that 

+ HVutlli* <c J g\u\\u t \ (|V Ut | + |Vn | 2 + | U ||V 2 U |) da: 

+ C J g\uf\S7u\\Vu t \&x + C J p|wt| 2 |Vu|dx + C J |Vd||Vdt||Vtit|da; 

=./1 + J 2 + J 3 + Ja- (3.26) 


By using (12.41) . (13.141) . Holder’s and Gagliardo-Nirenberg inequalities, it follows that 

Ji <C\\g?u\\ L 6\\g?u t \\l4g?u t \\l 6 {\\Vu t \\ L 2 + ||Vu||| 4 ) + C\\giu\\l 1 4g^ut\\l 2 \\g^u t \\l e \\V 2 u\\ L 2 
<C(l + ||Vu|| 2 2 )||^ Mt ||i 2 (||^ Mt || L2 + ||V Ut || i2 )5(||V Ut || i2 + ||V U || 2 2 + ||V W || L2 ||V 2 U |U2 + ||V 2 u|| i2 ) 


ll^tlll 2 + + ll£ ,2u tlll 2 ) + C'( 1 + I|Vm||1 2 )||V 2 w|| 


2 . 
L 2 ’ 


J 2 +J 3 < C\\g 2 u\\i s \\S/u\\ L 4^u t \\ L 4 + C\\g 2 u t \\l 6 \\g 2 u t \\l 2 \\S/u\\ L2 
<Jl|Vu t ||| 2 + C^(t)(l + Wg^utWl.) + C(l + ||V w || 2 i 2 )||V 2 w,|| 


2 . 
L 2 f 


Ja <C\\Vd\\ L 4 Vd t \\ L 4 \ 7 u t \\ L 2 < -||Vu t ||| 2 + C||Vd|| L 2 ||Vd|| ff i||Vd t || Z/ 2 ||Vd t || ff i 


<gl|Vu t || i2 


1 


4C*2 + 1 


||V 2 d,|| 2 2 +C<^(f)||Vd t || 2 2 , 


where the positive constant C 2 is defined in the following (13.281) and (13.311) . Substituting the estimates of 
Ji(i = 1,2, • • • ,4) into (13.261) . and then using (13.231) . we obtain 


^ll^ll ! 2 + livu t ||| 2 

<C$P(t)(l + ||g 2 Mt || 2 2 + ||Vd t || 2 2 ) + 


1 


4 C 2 + 1 


|V 2 d t ||i 2 +C(l + ||Vu|| 2 i 2 )||V 3 d|| 2 i2 . 


(3.27) 
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Now, differentiating (II. ll) 3 with respect to time variable t , and then multiplying the resulting equality by 
dt, we deduce that 

^||dt||l 2 + ||W t ||| 2 <C J \u t \\Vd\\d t \dx + C j |Vd t ||Vd||ci t |dx + C J |Vd| 2 |d t | 2 dx 
— J5 + Je + J7, 

where we have used the condition -§^d t = 0 on 8Br. Using Holder’s and Gagliardo-Nirenberg inequalities, 
and (13.141) . we have 

J 5 <\\u t x-%\\ L *\\\7dx%\\ L z\\d t \\ L * < (\\e^ut \\ L 2 + \\'Vut\\L 2 )\\Vdx%\\ L 2 \\dt\\l 2 \\d t \\] I i 

+ C<t>Pme>u t \\h + ||dt||l 2 + ||V<it 11^2), 

Je+Jr < i||Vd t || L 2 + C'||Vd|| 2 4 ||d t || 2 4 < i||Vd t || L 2 +C||Vd|| i 2||V 2 d|| i 2||d t || L2 ||d t || ffl 

<\\\Vdt\\ L *+C^{t)\\d t \\ 2 L ,. 

Hence 

l 2 + IIV'd*||5, 2 ^CallVutll^ + C'$' J (t)(||e 2 u t || 2 2 + ||d*||5,2 + IIVd*||5,2). (3.28) 

Differentiating (13.81) with respect to time variable t yields 

Vd tt - A= -V(u • Vd)t + V(|Vd| 2 d) t . (3.29) 

Multiplying (13.291) by Vd*, and integrating the resulting equality over Br, it follows that 

^||Vd i ||| 2 + ||V 2 d t || 2 2 + ^||Vd t || i 2 (aBji) <cj |Vuf ||Vd||Vd t |dx + C J |Vu||Vd t | 2 dx 
+ C J \u t \\V 2 d\\Vd t \dx + C J \Vd\ 2 \d t \\V 2 d t \dx + C J |Vd||Vd t ||V 2 d t |dx 
= Jg + Jg + J 10 + Jll + J\2- (3.30) 

By virtue of Holder’s and Gagliardo-Nirenberg inequalities, and (13.141) . we have 

Js <||Vu t || ia ||V*IU‘IIVd|U* < i||Vu t || 2 L 2 +C||Vd|U 2 ||Vd|| ff i||Vd t |U2||Vd t || ff i 
<^l|Vu t || 2 i2 + |||V 2 d*||2 2 + C^(t)\\Vd t \\ 2 L2 , 

J 9 <C'||Vu|| L 2||Vd t ||i4 < C , ||V«|| L 2||Vd t || ia ||Vd t || ff i < i||V 2 d t || L 2 + C^(t)\\Vd t \\ L 2, 

J 10 <c J \utX-^\\V 2 d\h$\V 2 d\^\Vd t \dx < C|| W4 x-t|| i4 ||V 2 dxt||J 2 ||V 2 d||| 2 ||Vd t || L 4 
<C||u t x-t|| L4 ||V 2 dxf||i 2 ||V 2 d||| 2 ||Vd t ||i 2 ||Vd t ||| 1 
<gl|V 2 dt|||2+G|| Ut x-t|| 2 4 + G||V 2 dxf|| 2 2 ||V 2 d|| 2 2 ||Vd t || 2 2 
<5l|V 2 dt|||2+G(||^^|| L 2 + ||Vu t || 2 2 ) + G||V 2 dxf|| 2 2 ||V 2 d|| 2 L2 ||Vd t ||| 2 , 

J 11 <|l|V 2 d*||i2 +C'||Vd||i 8 ||d t ||i 4 < i||V 2 d t ||i 2 + C||Vd||| 2 ||Vd||| 1 ||d t || i 2||d,|| ffl 
<|l|V 2 d f ||i2 +c^(t)(||d t ||| 2 + ||Vd t ||| 2 ), 

J 12 <^l|V 2 d t ||i 2 + C||Vd|| 2 4 ||Vd t || 2 4 < |||V 2 d*||i 2 + C^(t)\\Vd t \\ 2 L2 . 
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Inserting the estimates of Jt{i = 8 ,9, • • • ,12) into (13.301) . it yields that 
|l|Vci t ||| 2 + ||V 2 d t ||| 2 

<Ci{\\e^u t \\ L * + ||Vu t ||| 2 ) + C(&(t) + ||V 2 dit ||| 2 ||V 2 d||| 2 )(|| ( I t ||| 2 + ||W t || 2 2 ). (3.31) 

Finally, multiplying (13.271) by 2 C 2 + 1 and adding the resulting inequality with (13.281) and (13.311) . we have 

— (( 2 C 2 + l)||f ? 2 ut || 2 2 + WdtWni) + ||VMt||L 2 + -\\^d t \\ 2 H i 

<C(&{t) + ||V 2 ^f ||| 2 ||v 2 d|| 2 2 )(l + nautili. + HVdtl&O + (7(1 + ||Vu|| 2 2 )|| V 3 d||| 2 . 

Multiplying the above inequality by t, we obtain (13.241) after using Gronwall’s inequality and (13.101) . This 
completes the proof of Lemma 13.61 □ 


Lemma 3.7 Let (g,u,P,d) and T\ be as in Lemma 1,7.51 Then there exists a positive constant j3 > 1 such 
that for all t £ (0,Ti], 

sup (r|| V 2 m|| 2 2 + r|| V 3 d || 2 2 + r|| VP|| 2 2 ) < <7 exp j<7 exp j<7^ ^(r)drJJ . (3.32) 

Proof. Applying H k {k = 2,3) estimates of elliptic equations (see [13]), it is easy to deduce from (11.11b that 


||v 3 d || 2 2 <C||Vd t || 2 2 + C||M|V 2 d||| 2 2 + <7|||Vu||Vd||| 2 2 + C'llVdllle + C|||Vd||V 2 d||| 2 2 . 


By virtue of (13.141) . Holder’s and Gagliardo-Nirenberg inequalities, it follows that 

IIM|v 2 d||| 2 2 <C\\ux~HUV 2 dx%\\ L 4V 2 d\\ L . < c\\w^dxm +c\\ux-^\\t,\\v 2 d\\l< 

<q|v 2 d®t|^ + c -(||^ u || 2 3 + ||v u ||i 2 ) 2 || V 2 d|| ia ||v 2 d|| H i 

<^l|V 3 d||| 2 + C||V 2 da)f ||| 2 + (7(1 + ||Vu||| 2 )||V 2 d||i 2 , 

|||Vu||Vd|||| 2 <G||Vu|| 2 4||V 2 d|| 2 4 < ||V U || i2 ||V«|| ff i||Vd|U 2 ||Vd|| ir i 
<J||V 2 u||| 2 +C||Vu|| 2 2 (||Vd||l 2 + ||V 2 d||i 2 ), 

[|Vd||«e + |||Vd||V 2 d||| 2 2 < (7||Vd||| 2 1|Vd||^i + C||Vd|| 2 4||V 2 d|| 2 4 

<c\\Vd\\l 2 II Vdll^x + C7|| Vd|| i2 II Vdllffi II V 2 d|| i2 II V 2 d|| ff i 
<Jl|V 3 d||! 2 + G(l + ||Vd||l 2 + ||V 2 d||l 2 ). 

Hence 

l|V 3 d||| 2 <i||V 3 d|| 2 2 + J||V 2 u|| 2 2 + C||Vd t || 2 2 + C\\V 2 dx% II % + C{ 1 + ||Vu||i 2 2 + ||Vd||^), 
which together with (13.231) shows that 

||V 2 u||| 2 + ||VP||| 2 + ||V 3 d||| 2 

<C{\\e^u t f L , + llVdJi, + l|V 2 d.Tt|| 2 2 ) + C( 1 + || Vu||^ 2 + ||Vd||&)- 


Then, multiplying the above inequality by t, one obtains from (13.41) . (13.101) . (13.221) and (13.241) that 
sup (r||V 2 h || 2 2 + t||VP || 2 2 + r||V 3 d||^ 2 ) 

0 <T<t 


< (7 exp | (7 exp | (7 J ^(rjdr j | + C ^1 + exp |(7 J ' 
<C exp | C exp | C J ^(t)At 11, 


which completes the proof of (13.321) . 


□ 
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Lemma 3.8 Let (g,u,P,d) and T\ be as in Lemma \3.5\ Then there exists a positive constant /3 > 1 such 
that for all t £ [0, T±], 

sup ||£>5 a ||Li n ffi n wi>9 < exp |Cexp |Cexp |$ /3 (r)dr||| . (3.33) 

Proof. Inspired by [15] . we should first derive the following estimates 

I (l|V 2 n||Zf + IIVPH3 1 + r||V 2 u|| 2 g+r||VP|| 2 ,)dr<Cexp{cexp{c jfV(r)dr}}. (3.34) 
In fact, from (13.221) . Holder’s and Gagliadro-Nirenberg inequalities, we have for p = q, 


|V 2 zt|| L , + ||VP|| L9 <C\\ Q u t \\ Lq + C\\gu ■ Vu|| L5 + C\\ div(Vd © Vd)|| L * 

<G||^ t || L9 +C'||pu|| i2? ||Vu|| i 2, + C||Vd|| i 2,||V 2 (I|| i 2 9 


2(A-g) 


A(g—2) 


<\\gu t \\^\\gu t \\^ +C^{t)(l + \\V 2 u\^~ q +\\^d\\\^) 

2(A-g) A(g —2) l_i 1—1 

<C{\\ e ->u t \\W^ \\Vu t \\W^+ \\^u t \\ L 2) +C&m +l|V 2 u|| i 2 i + ||V 3 ( I|| i2 ^), 
where A > q will be chosen later. The above inequality together with (|3.10j) , (|3. 1 1|) and (|3.24 1) ensures that 


||V 2 u| 


L<? 


IIVPII 




dr 


pt . . (A g)(g+l) A(g_2)(g+1) /*£ q-j-i 

<C J (j\\g^ut\\ 2 L 2 j qHx ~ 2) (t||Vu t ||| 2 ) 2,2(a - 2) dr + CJ \\g^u T \\ L l At 


+ C ^(t)(l + ||V 2 u||^r + ||V 3 d||^r)dr 
Jo 

(A —g)(g+l) 


<C sup 
0 < 


Mg-Wg+l) 


up jr||p2 Mr ||| 2 l q {X 2> / t 2g (r||Vu r || 2 2 ) 2 ' j2(> '- 2) dr 

_T<t ^ ' JO 

+ c/ (4 >/3 (t) + ||0 2 u r || 2 2 + ||V 2 u||| 2 + ||V 3 d|||g)dr 

Jo 


<Cexp |Cexp 


{ C / V( 


r)dr 


i‘( 


g(g+i)(A — 2 ) 


1+ / IT 2g(A-2)-A(g-2)(g+l) + T11 Vw r || ^2 ) dr ) . 


By selecting A = 2<j 2 , we have 2 g(A ’^( ( g- 2 j( g + 1) = and then 


|V 2 M || L , + 


IIVPII L<? <Cexp |cexp J 4> /3 (r)di 

<C exp | C exp | C J $ /3 (r)dr|| . 


■7 J + q -g-i 


1+ / r + r||Vu T ||2, 2 J dr 


Similarly, 


(' r l|V 2 u|||g+T||VP|||g)da: < c j Q (^WQ^rWh)^ (r||Vu T ||| 2 )V-r d r (A = 2g 2 ) 

+ C [ T\\giut\\ 2 L 2 dT + C [ (^(r) + ||V 2 m||| 2 + ||V 3 d||| 2 )dr 
J 0 Jo 

<C J* (r||^||| 2 + T||VV|| 2 2 ) dr + cj\* f> (T) + ||V 2 u|| 2 2 + ||V 3 d|| 2 2 )dr 


<C exp | C exp | C J $ /3 (r)dr j> j> . 


Thus we obtain (13.341) . 
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Next, it is easy to derive from CCEDi that gx a satisfies 

d t {gx a ) + u ■ X7{gx a ) — agx a u • V lnx = 0, 

which implies for all r G [2 ,q\, 
d 


dt 


\\V(gx a )\\ L r <C( 1 + ||Vu||l- + ||u • Vlni||^)||V(px a )|| L . 

+ C'||^ o || i »(|||V U ||Vlnx||| i . + ||H|V 2 lnx||| i 0 

<C(&(t) + \\V 2 u\\L*nL*)\\V(gx a )\\ L r- + C\\gx a \\ L ~>(\\Vu\\ L r + \\ux~% \\ L *r\\x~i ||^) 


<C(* p (t) + ||V 2 M || i2ni9 )(l + ||V(e®°)|| L r + ||V( e ®“)||iO. 
where we have used (13.141) . (13.151) and the following calculatiorj 3 ] 


(3.35) 


II ux d \\ L °° <C{6){\\ux d \\ L *) + ||V(ux ^)||£ 3 ) 


<C{5){\\ux \\ L 4 + ||Vu||£3 + ||ux H^Hx Vx^^) 

<C(6)($P{t) + ||V 2 u|| L 2 ), for 0 < <5 < 1. 

Utilizing (13.341) . (13.111) and Gronwall’s inequality, it yields from (13.351) that 

sup ||V(px a )|| L 2 niP < exp jcexp jcexp |c J <f>' 9 (r)dr| 11 . 

This, along with (13.151) gives (13.331) . This completes the proof of Lemma [3781 
Now, we can deal with the proof of Proposition 13.31 

Proof of Proposition 13.31 It follows from (13.3l) . (13.41) . (13.51) . (13.101) . and (13.331) that 
$(£) < exp |cexp |cexp J $^(r)dr| ||, tG(0,Ti]. 

By using standard arguments yields that for M = e CeCe and T 0 = min{Ti, (CM' 3 ) -1 }, 

sup $(£) < M, 

0 <t<T 0 


(3.36) 


□ 


which together with (13.41) . (13.51) . (13.101) . (|3.11l) . (13.241) and (13.341) ensures (13.21) . This completes the proof of 
Proposition 13.31 □ 


3.2 The proof of Theorem 13.11 

With the a priori estimates obtained in Subsection 3.1, we shall give the proof of Theorem l3.ll Let {go, uo, do) 
be as in Theorem l3.1I Then it follows from (11.51) that there exists a positive constant Nq such that 


f g 0 {x)dx > 7 / g(x)dx = 7 . 
Jb„„ 4 J r 2 4 


1 Bn o “* 4r 2 

Now, we construct go = do + R~ l e~ l x l 2 , where 0 < go £ Cq°(R 2 ) satisfies 

fe No do(x)dx > 

gffx a —> g 0 x a in L 1 (R 2 ) n H 1 {M. 2 ) D fP 1 , 9 (M 2 ), as R —> 00 . 


!We notice that the inequality (13.361) still holds if we replace the integral domain Br by M 2 . 
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Noticing that do G T 2 (R 2 ;§ 2 ), SJdoX^ G T 2 (R 2 ) and V 2 do G T 2 (R 2 ), standard arguments (cf. [29] and 
also [13]) yield that there exists d R G -ff 3 (R 2 ;§ 2 ) such that d R = no outside Br for some constant vector 
no G § 2 and 

lim \\X7d R - Vd 0 ||if 2 (R 2 ) = 0, lirn \\Vd R x% - Vd 0 x% || L 2 ( R 2 > = 0. 

Ft—too R—too 

Since Vrto G T 2 (R 2 ), choosing v R G Cq°(-Br)(* = 1,2) such that for * = 1,2, 

lim ||«f - 3i«o|U 2 (R 2 ) = 0 . 

-R->oo 

Let u R be the unique smooth solution of the following elliptic problem : 

-A u R + g R u R + VP R = ^$h R - d iV f, 

div u R = 0 , 

u R = 0 , 

where h R = (pgUo) * jh with js being the standard mollifying kernel of width 6. Extending u R to R 2 by 
defining 0 outside Br and denoting it by u R , then by the same arguments as those for the proof of Theorem 
1.1 in [T9] (see also [35]), we obtain that 

“ u o)IU 2 (R 2 ) + \\(So)^o - 0o u o||l 2 (r 2 )) = 0. 

Then, by virtue of Lemma 12.11 the initial-boundary value problem (11.11) and (12.21) with the initial datum 
(g R ,u R ,d R ) has a classical solution (g R ,u R , P R ,d R ) on Br x [0, Tr]. Furthermore, Proposition 13.31 shows 
that there exists a To independent of R such that ([3.21) holds for ( g R , u R , P R , d R ). Extending (g R , u R , P R , d R ) 
by zero on R 2 \Br and denoting it by 

(g R ±v R g R ,Z R ,P R ,d R ), 


in Br, 
in Br, 
in 8B r , 


with ipR as in (13.121) . it deduces from (13.21) that 


sup ^||(p' R ) 2 W' R || i 2 ( R 2 ) + ||V*t fl || i 2 ( R 2 ) + ||V 2 d|| i 2 ( R 2 ) + ||Vdx 2 || L 2 ( R 2 )^ 


0<t<T 0 


< 


sup (||( 0 fl )^|U 2(Bfl) + \\Vu r \\l HB r) + ||V 2 d|U 2 (Bfl) + IlVdsf |U 2 (Bb) ) < C, (3.37) 


0 <t<T 0 


and 


sup ||p a; 0 ||L 1 (R 2 )ni<' x ’(R 2 ) < C sup ||e x a \\ L i {BR)nL ^ {BR) < C. 
o <t<T 0 o <t<T 0 


(3.38) 


Similarly, it follows from (13.21) that for q > 2, 

sup t 2 f Ke^) 2 u R ||l 2 ( R 2 ) + ||df ||_f*i(R 2 ) + || V 2 rr R || i 2 ( R 2 ) + || V 3 f7 ; ' | B 2 ( R2 ' 
0<t<T o K 


rTo 


+ 

+ 


(||(e fl ) 2 uf ||| 2 ( R 2 ) + ||V 2 M y 'j|| 2 ( R 2 ) + ||df || B i( R 2 ) + ||V 3 d fl || B 2 ( H 2 ) + ||V“d fl a ; 2 ||| 2 ( R 2 )^ dt 


r T o / 9+i _ 

ll^ 2 ^!! li{ r 2 ) + i ll^ 7 “ u " R |ll<i(R 2 ) + i||Vnf||| 2 (R 2 ) + i||Vdf|| B i( R2 ) j < C. 


(3.39) 


Next, for p G [2, q], it follows from (13.21) and (13.331) that 


sup ||V(e fi a; a )|| B p( R 2 ) <C sup (\\V(g R x a )\\ L p {BR ) + R 1 \\g R x a \\ L p {BR )) 

0 <t<T 0 0 <t<T 0 
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(3.40) 


<C sup \\g R x a \\ H i( BR ' >nw i,p( BR ) < C, 

0 <t<T 0 


which together with (13.361) and (13.21) ensures 

rTo 


\\Q?x a \\ 2 L P(R 2 )dt < C [ \\\u R \\Vg R \x a \\ 2 LP(BR) dt 

Jo 


<c ( T ° ||^ 1 -l|o 0 (BR) ||v^ a ||| P(BR) dt < c. 


(3.41) 


JO 

By virtue of the same arguments as those of (13.231) and (I3.34|) . it follows that for q > 2 

i ~ r T ° ( ~ ~ S±i \ 

sup V s ||VP' R || B 2 (R 2 ) + / II VP h |||, 2 (R 2 ) + ||VP*y (R2) < C. (3.42) 

With all these estimates (13.371) (13.421) at hand, we find that the sequence ( g R ,u R ,P R ,d R ) converges, up 
to the extraction of subsequences, to some limit (p, it, P, d) in the obvious weak sense, i.e., as R —>• oo, we have 


g R x —»• gx in C(Bn x [0,To]) for all N > 0, 

g R x a ->• gx a weakly * in T°°(0, T 0 ; H^R 2 ) n W^^R 2 )), 

Vd R x% ->• Vdx% weakly * in L°°(0, T 0 ; T 2 (R 2 )), 

(g R )hu R ->■ 0 5it, Vd R ->■ Vu,V 2 d R V 2 d weakly * in T°°(0, T 0 ; T 2 (R 2 )), 

V 2 y« V 2 u, VP fl ->• VP weakly in (0, T 0 ; T?(R 2 )) H P 2 (R 2 x (0,T„)), 
df —*■ dt weakly in T 2 (0,T o ; P^R 2 )), 

< V 2 d%^ ->■ V 2 dx% ,S7 3 d R V 3 d weakly in T 2 (M 2 x (0,T o )), (3.43) 

t?V 2 u R -+t?V 2 u weakly in T 2 (0, T 0 ; T«(R 2 )), weakly * in T°°(0, T 0 ; T 2 (R 2 )), 

£{q r )^u? g^Ut,t^\7P R UX7P weakly * in T°°(0, T 0 ; T 2 (R 2 )), 

f5df -a t \dt weakly * in T°°(0, T 0 ; P^R 2 )), 

t*V 3 d R t?V 3 d weakly in T°°(0, T 0 ; T 2 (R 2 )), 

t^\7u R —- tiVu* weakly in P 2 (R 2 x (0,To)), 

f^Vdf-^Vdt weakly in P 2 (0, To; P 1 (R 2 )), 

with 

gx a G T°°(0,T o ; T 1 (R 2 )), inf f g(x,t)dx>^-. (3-44) 

o<i<T 0 ,/ B2JVo 4 

Then, letting P —> oo, standard arguments together with (13.431) and (13.441) yield that ( g,u,P,d ) is a strong 
solution of system ED ED on R 2 x (0, To] satisfying (11.81) and (11.91) . Indeed, the existence of a pressure P 
follows immediately from the equations ED 2 an d ED 4 by a classical consideration. Thus we complete the 
proof of the existence part of Theorem 13.11 

In what follows, we shall prove the uniqueness of the strong solutions. Let (gi, u\, Pl, di) and (q 2 , U 2 , P 2 , d 2 ) 
be two strong solutions satisfying ED and ED with the same initial data. Let g = gi — 02 , u = u\ — 112 , 
P = Pi — P 2 and d = di — d 2 , it follows that 

gt + U 2 ■ Vp + u • Vpi =0, 

g\Ut + gi Hi • Vu —Au = giu ■ Vu 2 ~ g(u 2 t + U 2 ■ V 119 ) — VP — div(Vd © Vdi)— div(Vd 2 © Vd), 

< _ _ (3.45) 

dt — Ad = — U 2 ■ Vd — u • Vdi + |Vd2| 2 d + V(di + d2)Vddi, 

div u = 0, 
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for ( x , t) G K 2 x (0, T 0 ] with 

g(x, 0) = u(x, 0) = d(x , 0) = 0. 

Multiplying (13.451) x by 2 gx 2r for r G (l,a) with a = min{2,a}, and then integrating by parts over R 2 yield 


— I \gx r \-&x<C\\u 2 xi || L =o (R 2 ) ||ea: r ||£ 2(R2) + C'||£>a; r ||L 2 (R 2 ) ||Ma; 


--(a-r)i 


I 2 q |foia; a || 2q 

i (,-2)(a-r) ( R 2) i5 -(,-2)(5-r) (R2) 


<C{1 + ||w2||vv 1 >9( R 2))||pa: , ’|| L 2( R 2) + C||£ | a; r '||L2( R 2)(||Vu|| i 2( R 2) + \\qI w|k 2 (R 2 ))) 


due to Sobolev’s inequality, m, Lemma [2.31 and (13.361) . This combined with Gronwall’s inequality shows 
that for all 0 < t < T 0 , 


(3.46) 


II^1l2 (R 2) <C (||Vm|| L 2 (R 2 ) + ||£< 1 2 M||L 2 (R 2 ) )dr. 


Multiplying (I3.451) 3 with d , and then integrating by parts over R 2 , it follows from (11.81) and Lemma [2731 


that 


~ [ \d\ 2 dx+ f \Vd\ 2 dx 

z at J R 2 JR 2 

<[ \u\\S 7 d l \\d\dx+ f \S 7 d 2 \ 2 \d\ 2 dx+ ( |V(di + d 2 )||Vdj|d|da; 

J R 2 J R 2 Vr 2 

<C\\UX 4 ||_L 4 (R 2 ) WVdlX 2 |||,2( R 2)||V(il||| ;2 ( R2) ||dj|L4( R 2) + (7||Vd2|| 2 4( R 2)||dj||4( R 2) 
+ C'dl Vc?l ||i4( R 2) + ||Vd2||L 4 ( R 2))||d||L 4 ( R 2)||Vd|| L 2( R 2) 

<e:|| Vrf||^2( R 2) + e\\ux 4 ||| 4 ( R 2 ) + ( 7 ||Vdi 5 2 |||2( R 2) || Vdi||| 2 ( R 2)||d|||2( R 2) 

+ C'(e)(l|Vdi||i2 (R 2 ) ||V 2 di||| 2(R 2 ) + ||Vd2||| 2 ( R 2) ||V 2 d 2 ||l2( R 2))||<i|||2( R 2) 

<£r|| + £(||0 iu|Il 2 (r 2 ) + ||Vm||| 2(r2) ) + C'(£)||d||| 2 ( R 2), 


(3.47) 


where we have used |di| = 1 , and the fact that the divergence free condition divu .2 = 0 ensures the identity 
f R2 u 2 ■ Vd • ddx = 0 in the first inequality. 

Next, multiplying (13.451) .-, and (13.45IL by u and —A d respectively, and adding the resulting equations 
together, we obtain after integration by parts over R 2 that 


1 d 

2 dt 


( gi \u\ 2 + \Vd\ 2 )dx + / (|W | 2 + |V 2 d| 2 )dz 
J J R2 

<C'||Vu 2 ||l~(r 2 ) / (gi\u\ 2 +\Vd\ 2 )dx + | 0 ||w|(|u 2 t| + |w 2 ||Vu 2 |)d;z; 


|w||Vdi||Ad|da; 


|Vd 2 | 2 |d||Ad|da; ■ 


(|Vdi| + |Vd 2 |)|Vd||Vu|dx 

J R 2 

(|Vdi| + |Vd 2 |)|Vdj|Adjdz 


—(711 Vm 2 ||loo( R 2 ) / (gi|tt|“ + |Vd| 2 )d.T +//i + J / 2 + • • • + 1 / 5 . 

Jr 2 

By using Holder’s inequality, Lemma [2731 and (13.461) . we have for r G (l,a), 

Hi <C\\gx r \\ L 2 ^^\\ux 5 ||l 4 ( R 2 ) (||u 2 tx 5 ||l 4 (r 2 ) + ||Vu 2 ||l~(r 2 )||m25 2 ||i 4 (R 2 )) 


(3.48) 


<£ 


\ UX 2 lli 4 (R 2 ) + C( e ) (\\S2 u 2t|| 2 2( R 2) + ||Vu 2 t|| 2 2( R 2) + IIV ^2 II 


L°°(R 2 ) j II Q x ' lli 2 (R2) 

t 


— £ ( 11 01 U llL 2 (R 2 ) + ll^ 7M llz, 2 ( R 2 )) + C'(£)(l+t|| Vu 2 t|lL 2 ( R 2 )+ t||V W2||l9( R 2)) / (|| 01 n||i 2 ( R 2 )+ ||Vu|| L 2 ( R 2 ))dr. 


'0 


For II 2 and II 5 , we derive from Holder’s and Gagliardo-Nirenberg inequalities, and (11.81) that 
II 2 + II 5 <C'(||Vdl|| i 4( R 2) + ||Vd2||L 4 (R2))||Vd|| L 4 (R 2 ) (||Vu|| i 2 (R 2 ) + ||V 2 d|| i 2 (R 2 ) ) 
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< £ (II Vu IIl 2 (r 2) + liv-dllla^)) + C'(e)||Vd||i2 (Ra) . 

By using Holder’s and Gagliardo-Nirenberg inequalities, and Lemma \2 . 31 again, we can estimate II 3 as 

II 3 <C\\UX ~2 ||i-l(R 2 )||VdlX 2 ||i4( R 2) || V 2 d|| L 2 ( R 2 ) 

< £ l|V 2 d||l2 (R 2 ) + C(e)( 1 + ||V“dix 2 || 2 2( R 2))(||£>i m||| 2 (R 2) + ||Vzt||| 2 ( R 2)). 

The term II 4 can be estimated as follows 

Hi <C\\Vd 2 ||£ 8 (R 2 ) |M||l4( R 2) || V 2 d|| L 2 (R 2 ) < e||V 2 dj| 2 2 ( R 2 ) + C(ff)||d||^-l( R 2 ), 

owing to Holder’s and Gagliardo-Nirenberg inequalities, and (11.811 . Inserting these estimates //,;(* = 1, 2,••■ , 5) 
into (13.4811 . adding the resulting inequality with (13. 4711 . and then choosing e suitably small lead to 

G'(t) < <7(1 + ||Vw 2 ||l”(r 2 ) + ||V 2 dia ;2 ||| 2 (R 2 ) + t\\Vu 2 t\\ 2 L 2 ^ +t\\S/ 2 U 2 \\ 2 LQ ^ M 2 ^)G(t) (3.49) 

where G(t) is defined as 

G(t) = ||^ti|||2( H 2) + ||Vd||^-i( R 2) + f (|| Vw|| 2 2 (R 2) + || V^d|| 2 2( R2 ) + ||<?i m|| 2 2 )dr. 

Jo 

The inequality (13.491) together with Gronwall’s inequality and (11.81) implies that G(t) = 0 for all 0 < t < T 0 . 
Hence, we have u = 0 and d = 0 for almost everywhere (x, f) £ R 2 x (0,T o |. Finally, one can deduce from 
(13.4611 that g = 0 for almost everywhere (x,t) el 2 x (0,To]. This completes the proof of Theorem 13. II □ 

4 Proof of Theorem 11.2 

Before going to do it, let us recall the following rigidity theorem, which was recently established in [16]. 

Theorem 4.1 Let Eq > 0 as 0 , and Mo > 0. There exists a positive constant uj = w(eo,Mo) £ (0,1) such 
that the following holds: 

If d : R 2 —S 1 , Vd £ H 1 ) K 2 ) with || VcZ||x ,2 < M 0 and d 2 > £o, then 

ll^^lli 4 (R 2 ) < (1 — d))|| Ad||| 2 ( R 2 ). 

Consequently for such maps the associated harmonic energy is coercive, i.e., 

||Ad+ |Vd| 2 d||| 2 ( R 2 ) > — (||Ad||| 2 ( R 2 ) + ||Vd||| 4 ( R 2 )) . (4.1) 

In what follows, for p G [1, oo] and k > 0, we denote 

f fdx = [ fdx, LP = L P {M. 2 ), W hp = LF 1 ,P (M 2 ), H k = W k ’ 2 , 

J J R 2 

and denote by || ■ ||x the norm of the X(R 2 )-functions, for simplicity. 

4.1 Lower order estimates 

In this subsection, we shall establish some necessary a priori lower order estimates for strong solutions 
( g,u,P,d ) to the Cauchy problem of system (11.11) (11.21) . Let T > 0 be a fixed time and ( g,u,P,d ) be the 
strong solution to system (11.11) (11.21) on E 2 x (0, J 1 ] with initial data {go, uq, do) satisfying (11.31) and (11.51) (11.71) . 
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Lemma 4.2 There exists a positive constant C depending only on £o, UpolUini 00 ) Ildo^oIlL 2 ; ||Vuo||l 2 , 
IIVc^oHz . 2 and ||V 2 do||z , 2 such that 


SU P (IMIzPnz , 00 + l|Vw|| 2 2 + || V“d|| 2 2 ) + f (lie 2 w||i ,2 + ||V 3 d|| 2 2 ) dt < C, (4.2) 

o <t<T Jo y 7 

where u = Ut + u ■ Vu. Furthermore, we have 

sup t (||Vu ||| 2 + ||V 2 d||| 2 ) + [ t (||e ^||| 2 + ||V 3 d||£ 2 ) dt < C. (4.3) 

0 <t<T Jo K 7 

Proof. We first notice that, the basic energy law (13.31) and the identity (13.61) still holds in R 2 . Under the 
assumption & then it follows from the standard maximum principle method that d 3 > e 0 (see |31|). Hence, 
C3D and (13.31) together with (14.11) in Theorem 14.11 it follows that for all t > 0 

sup (||e 2 u ||| 2 + ||Vd||| 2 ) + - /" (||Vw || 2 2 + ||Ad || 2 2 + ||Vd||^ 4 ) dr < ||pq M 0 III 2 + ||Vc?o||| 2 . (4.4) 

0 <T<t v 7 ^ Jo 

Now, multiplying 0> 2 by u , and integrating over R 2 , it follows that 


/ 


e|ft| 2 da; 


Aw ■ udx — 


VP • wda; — 


div(Vd © Vd) • iidx 


A/ 1+ / 2 + J 3 . 


(4.5) 


By using the definition of u , and the Gagliardo-Nirenberg inequality, we have 


1 1 = J A u ■ (ut + u ■ Vw)da; = — J Vu • (Vut + V(w • Vw))da; 
2 . r 

di Uj diUkdkUjdx 


5sll v “Hb- £ 


2 dt 
1 d 


i,j,k=l' 


1 d 


< ~^W\h + H v “ll£» ^ ~2 ^U Vu ll£ 2 + C||Vw||| 2 ||V 2 w|U 2 , 


where we have used the divergence free condition <hhd 4 in the first inequality. By integrating by parts, and 
then using the divergence free condition (11. ID t and the duality of P 1 (R 2 ) and BMO( R 2 ) (see p25], Charpter 
IV), it follows that 


1 2 


VP ■ (ut + u ■ Vu)dx 


PV{u ■ Vw)dx 



2 

PdiUjdjUidx < ||P||sMol|3iU • Vui\\ H 1 . 

i= 1 


Notice that div(c*iw) = d, divw = 0 and V 1 - • (Vui) = 0 for i = 1, 2. Hence, the above inequality together with 
Lemma 12.61 ensures that 

2 

I 2 < C^2\\P\\BMo\\diU ■ Vui\\ n i < CIlVPIl^llVwll 2 .. 

2 = 1 


To bound the term / 3 , integrating by parts together with (13.81) . Holder’s and Gagliardo-Nirenberg inequalities 
gives 


13 = J (Vd © Vd) • Viqda; — J div(Vd © Vd) • (u ■ Vw)dx 

=—J(Vd © Vd) • Vwdx — J(Vdt © Vd) • Vudx —J(Vd © Vdt) • Vwd.T —J div(Vd © Vd) -(u • Vw)da 
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= ltj 0 ' Vudx ~f K AW " V 0 ' Vd ) + V (l Vd l 2fi )) © Vrf ] • Vudz 

- J[VdQ (AVd - V(u • Vd) + V(|Vd| 2 d))] ■ S7udx -j div(Vd © Vd) ■ (u ■ Vu)da; 

2 

=|/(Vd ©Vd) -Vudz -J[(AS7d - Vu-Vd +V(|Vd| 2 d)) ©Vd] -Vudx+ E / Ukdkdid ■ djddjUidx 

i,j,k=l 

2 

- /[Vd © (AVd - Vu•Vd + V(|Vd| 2 d))] • Vudz + E Ukdid ■ dkdjddjUidx 

J i,j,k=l ^ 

+ did ■ djddj(ukdkUi)dx 

i i j 7 k=l 

=j t J (Vd © Vd) • Vudx -J [(AVd - Vu • Vd + V(|Vd| 2 d)) ©Vd] • Vuda 

2 

- /[Vd 0 (AVd-V.-W + V(|Vd| 2 d))] ■ V»d* +E/W- 


i,j,k= 1 1 

<1 |(Vd © Vd) • Vudx + C||V 3 d|| L 21| Vd|| L aIIVu|| L 3 + C||Vd|| 2 a IIVu||£ 3 
+ CllVdlliallVulUa + C'||V 2 d|| i 3||Vd|||a||Vu|| I/ 3 

d 
dt 
d 
dt 
d 
' dt 


<— j (Vd © Vd) • Vudx + -|| V 3 d ||/ 2 + C\\ Vu ||£ 3 + CfVd||£ s + C'lrV 2 d||/ 3 1| Vd||£a (with e > 0) 


< 


<- 


lie 


J (Vd © Vd) • Vudx + |||V 3 d||| 2 + C||Vu||£ s + C'lrVd||) 

/(Vd © Vd) • Vudx + |||V 3 d||| 2 + C , ||Vu||£ 2 ||V 2 u || i 2 + C'||Vd|| 2 2 ||V 2 d||i 2 . 

Inserting the estimates of /,; (i = 1,2,3) into (14.51) . and then using (14.41) . it follows that 

J t Q||Vw||| 2 - |(Vd © Vd) • Vudx^j + \\ e ^u \\\ 2 
<fl|V 3 d||| 2 +C||V 2 d||i 2 +C(||V 2 U || i2 + ||VP|| i2 )||Vu|| 2 2 . 

Notice that (g, u, P , d) satisfies the following Stokes system 

—Am + VP = — gii — div(Vd © Vd), i£l 2 , 
divw = 0 , i£l 2 

u(x) —> 0 , |ar| —>• oo. 

Applying the standard L p -estimate to the above system (see [ID]) ensures that for any p £ (1, oo) 

IIV 2 m||lp + IIVPHlp < C\\gu\\ LP + C|||Vd||Ad||| L P < C\\giu\\ LP + C|||Vd||V 2 d||| L P, 


(4.6) 


(4.7) 


where we have used the identity div(Vd • Vd) = Vd • Ad and (13.61) . Combining (14.61) and (14.71) together, it 
follows that 


dt 


B(t) + \\e^\\b <d|V d d || 2 i2 +C||V 2 d||l2 +C'||V «||£ 2 


|^u || 2 2 + |||Vd||V 2 d|||| 2 ) 


where 


<|l|v 3 d ||| 2 + C'||V 2 d||l 2 + C||V«||i a + £ -(\\g^ii\\h + ||Vd|| i 2 ||V 2 d|| 2 2 ||V 3 d||| i2 ) 
<e|| V 3 d ||| 2 + s\\giu\\l 2 + C\\S7 2 d\\\ 2 + C|| Vu|| 4 i2 , 

B(t) = i||Vtt ||£ 3 - J (Vd • Vd) • Vudx, 


(4.8) 
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satisfies 


j\\Vu \\ 2 l2 - C^dWl, < B(t ) < C||Vu || 2 2 + C||V 2 d || 2 2 
owing to Gagliardo-Nirenberg inequality, (TOl) and the following estimate 

/(Vd ■ Vd) ■ Vudx <|llV«||i a + C\\Vd\\ 4 Li < \\\Vu\\ 2 L2 + C||Vd|| 2 2 ||V 2 d || 2 2 

<^l|Vu |||2 + q|V 2 ri|| 2 2 . (by dm)) 

Now, multiplying da by —VAd and then integrating by parts over R 2 , it follows from Holder’s and 
Gagliardo-Nirenberg inequalities, 03) and 03) that 


lj f \\V 2 d \\ 2 L 2 + ||V 3 d ||| 2 = J V(« • Vd) • VAddx - J V(|Vd| 2 d)VAddx 

2 

= J (Vtt • Vd) • VAddx + ^ J mdidjd ■ didlddx — J V(|Vd| 2 d)VAddx 


2 


= J (Vu • Vd) • VAddx — ^ J dkUididjd ■ didkddx — J V(|Vd| 2 d) • VAddx 


i,j,k=l 


<G(|| V 3 d|| i2 || Vu\\ L3 ||Vd|| i6 + || Vu\\ L s || V 2 d|| 2 3 + || V 3 d|| i2 1| Vd|| 3 e + || V 3 d|| i2 1| V 2 d|| i3 1| Vd|| i6 ) 

<\ l|V 3 d||| 2 + C||Vu|| 3 3 + C||Vd|||a + C||V 2 d|| 3 3 + C||V 2 d|| 2 3||Vd|| 2 a 

<f l|V 3 d||| 2 + G||Vu|| 3 3 + C1|Vd||£a + C||V 2 d||i 2 

<|l|V 3 d||| 2 + G||Vu||| 2 ||V 2 u||i2 + C||Vd|| 2 2 ||V 2 d||l 2 + C'||V 2 d||l 2 

<f l|V 3 d||| 2 + |(||^u|| 2 2 + |||Vd||V 2 d||| 2 2 ) + C||V 2 d||i 2 + C'||V«||£ a 

<e||V 3 d|| 2 2 + e\\ e iu\\l 2 + C'||V 2 d||l 2 + C||Vu||i 2 . (4.9) 


Multiplying (14.91) by (Cl + 1), then adding the resulting inequality with (14.81) and choosing e suitably small, 
we obtain 

— (B(t) + (Ci + l)||V 2 d||| 2 ) + ||e a «||| 2 + ||V 3 d ||| 2 

<C(||V 2 d||4 2 + ||V«Hi a ) < C(||V 2 d||i 2 + ||V W || 2 i2 )(H(t) + (Ci + l)||v 2 d||| 2 ). (4.10) 

The above estimate (14.101) combined with energy inequality (14.41) . the definition of B(t), and Gronwall’s 
inequality ensures 

sup (||Vw || 2 2 + ||V 2 d||i 2 ) + f T (\\Q?u\\h + ||V 3 d|| 2 2 ) d t < C. 

0 <t<T Jo K ' 

This together with (13.61) ensures (14.21) . 

Finally, multiplying (14.101) by t, and then applying Gronwall’s inequality to the resulting inequality, it 
follows from (14.41) and the definition of B(t) yields (14.31) . This completes the proof of Lemma T4. 21 □ 

i 

Lemma 4.3 There exists a positive constant C depending only on £q, HgolliPnL 00 ) lldo^olli 2 ) ||Vuo||l 2 ; 
||Vdollx, 2 and ||V 2 do||z ,2 such that for i = 1,2 

sup f (\\giu\\ 2 L2 + |||Vd||V 2 d||| 2 2 ) + [ T (||Vd || 2 2 + |||Vd||AVd||| 2 2 ) dt < C, (4.11) 

0 <t<T v 7 Jo 

and 


sup 

o <t<T 


s (I|V 2 U || 2 2 + ||VP|| 2 2 ) < c. 


(4.12) 
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Proof. Operating d t + u • V to (ll.ll) 2j (j = 1, 2), one obtains by some simple calculations that 


dt(giij) + div(guuj) — Aiij 

2 

= ^ ^ — di(diu • Vuj) — div(diudiUj) — d t di{did ■ djd) — u ■ Vd^did ■ djd — d t djP — u ■ VdjP. 

i =1 

Multiplying the above equality by Uj, and then integrating by parts over R 2 , it follows that 

- — Hollis + ||V?i |||2 = - ^2 [di(di u • V Uj ) + di v(diudiUj)\ Ujdx 

i,j =1 ^ 

-%h a ‘ p+ »- - t(h md ■ w^ + /“ ■ va(M ■ 

J = 1 2 ,(7 = 1 X 


—J i + Ji + J 3 . 


(4.13) 


In what follows, we shall estimate each term on the right-hand side of (14.1311 term by term. We first notice 
that exactly the same arguments as Lemma 3.3 of paper 1361 . one has 

2 

h + J2<f t I ]T pa j i t< a i uid a: + c'(||P||t. + ||v«||i.) + i||VTi||| a . 

i,j =1 

For the term J 3 , by using d 3 > (HHD 4 and Gagliardo-Nirenberg inequality, it follows that 

J 3 = J diiijdidt • djddx + J diiijdid ■ djdtdx — J Ukdkdi(d,d ■ djd)ii 3 dx 


i,j =1 
2 


*d=l ’ 


i,j,k=l * 


= ^ f diiijdjd -di(Ad—u • Vd+|Vd| 2 d)d:r-|- ^ f diiijdid -dj{Ad—u • Vd+|Vd| 2 d)d:r 
i,j=l^ i ,j=V* 

— / Ukdkdi(did ■ djd)iijdx 

i,j,k,e= 1 J 

= E [diiijdjd ■ (Ad t d - dju ■ Vd + di(\Vd\ 2 d))dx - [ diiijdjd • dkdidukdx 

i,j=1 i,j,k,£= 1 

+ ^ / diiijdid ■ (Adjd - d 3 u ■ Vd + d 3 (\Vd\ 2 d))dx — E diiijdid ■ dkdjdukdx 

i,j=± i,j,k= 1 ^ 

+ y] / d l Ukdk(d l d ■ d 3 d)u 3 dx + ^ / ukdk(did ■ djd)diiijdx 

A A 7 ,_ 1 ^ 7 _ 1 J 


i,j,k =1 ' 




= ^ J diiijdjd • (Adid — diU • Vd + <9j(|Vd| 2 d))dx 
+ ^ / diiijdid ■ (Adjd — d 3 u ■ Vd + 5, (|Vd| 2 d))da; + E diUkdid ■ djddkiijdx 

A A 1 ^ A A 7 „ n 1 


*|J = 1 


i,j,k,Z =1 ’ 


<C'||Vu|| i 2 |||Vd||AVd||| i 2 +C||Vu|| L 2 ||Vu|| i 4 |||Vd | 2 || i 4 
— ^l|Vd || 2 2 + C|||Vd||AVd ||| 2 2 + C||Vii ||^ 4 + C|||Vd| 2 ||^ 4 . 

Inserting the estimates of Ji(i = 1, 2, 3) into (14.131) . one obtains 
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<j t I E ^ j « i a iUj d a; + C(||P||| 4 + ||Vn|||4 + |||Vd| 2 ||l4) + C|||Vrf||AVd|||| 2 

*,J=1 


(4.14) 


Now, inspired by the papers [34l(35ll37] . for oi, 02 € {—1, 0,1}, let us denote 

(Vd)(ai, 02 ) = ai<9id + 02^2 d, (Vu)(ai, 02 ) = ai<9iit + 0 , 28211 , u(a 1 , 02 ) = aiiti + a 2 U 2 , 

then it is easy to deduce from (13.81) that 

Vd t - AVd = -Vu • Vd-u ■ \7S7d + | Vd| 2 Vd + 2(Vd ■ d)VVd. 

it follows that 


Multiplying the above equality by 4VdA|Vd| 2 , and then integrating by parts over 

d_ 

dt 


: II V| Vo?| 21,2 


L 2 


2|| A|Vd| 2 " 2 


L 2 


= - 4y (Vu • Vd) • VdA|Vd| 2 da; - 4 j (u ■ VVd) • VdA|Vd| 2 dx 
+ 4,J |Vd| 2 |Vd| 2 A|Vd| 2 d.T + 8 ^Vd • d)VVd • VdA|Vd| 2 dz 
<C / |Vu||Vd||Vd||A|Vd| 2 |d:r — 2 /(it ■ V|Vd| 2 )A|Vd| 2 d:r 


+ C J Vd| 2 |Vd| 2 |A|Vd| 2 |d:r + C J |Vd|V|Vd| 2 |A|Vd| 2 |da; 

<C||Vu|| L 4|||Vd| 2 || i 4||A|Vd| 2 || i2 + C||Vu|| i 4||V|Vd| 2 || 2 f 

+ C|||Vd| 2 ||i4||A|Vd| 2 || i2 +C'||A|Vd| 2 |U 2 ||V|W| 2 || if ||Vd|U8 


<-||A|Vd| 2 ||i 2 + C(|||Vd| 2 ||£ 4 + ||Vit|| 4 4 ) + C||V|Vd|\ 

2 L3 


:| 2 p 

m2 II2 1 n(\\ I V7^7l 2 II4 1 ||Y7„.I|4 \ 1 n IIIV7^|2||| 

' 12 112 , /^/111 \~7 7 12 11 4 , ||V7„.||4 \ , ^Ilv72j||4 


■C'llVdllla 

r|2||| 


<-||A|Vd| 2 ||i 2 + C(|||Vd| 2 p4 + ||Vu||^) + C'|||Vdn|| 4 ||A|Vd| 2 ||| 2 + C||Vd|| 4 i2 ||V 2 d||£ 4 


<l|A|Vd| 2 ||£ 2 + C(\\\Vd\Tn + ||Vu||£ 4 ) + C'||V 2 d|| 4 i 4, 


(4.15) 


where we have used CQD 4 Holder’s and Gagliardo-Nirenberg inequalities, and (14.41) in the above estimates. 
Noticing that 

|||AVd||Vd|||| 2 <C||V 2 d||| j 4 + ||A|Vd(l,0)| 2 ||| 2 + ||A|Vd(0, l)| 2 ||| 2 + ||A|Vd(l, l)| 2 ||| 2 + ||A|Vd(l, —1)| 2 ||| 2 , 


and 


ll|V 2 d||Vd||| 2 2 < G(t) < C|||V 2 d||Vd||| 2 2 (4.16) 

with 

G(t) A ||V|Vd(l,0)| 2 || 2 2 + ||V|Vd(0,l)| 2 || 2 L2 + ||V|W(l, 1)| 2 1|| 2 + ||V|Vd(l,-1)| 2 || 2 2 . 

Thus, it follows from (14.151) multiplied by (C 2 + 1) that 

j t m + 1 )G(t)) + (C 2 + l)|||AVd||Vd||| 2 2 < C||Vu|| 4 4 + q|V 2 d|| 4 4 + C|||Vd| 2 p4, 
which combined with (14.141) ensures that 

j t m + |l|VA||i 2 + |||AVd||Vd|||| 2 < C\\P\\U + C||Vu||| 4 + C||V 2 d||i 4 + C\\ |Vd| 2 || 4 4 , (4.17) 

where 

1 r ^ 

F(t) A -\\giu \\ 2 L 2 + (C 2 + 1 )G{t) - J E PdjUidiUjdx 

»,j=i 


26 






satisfies 


\\\Q h n\\h + ~^Y~G(t) - CHVulli. < F{t) < \\ghi \\ 2 L2 + CG(t ) + C\\V U \\l* 


(4.18) 


owing to the following estimate 


2 

< ^ C\\P\\BMo\\diU • VuiUw 1 < C'||VP|| Z( 2 ||Vu |||2 (by Lemma HU) 

2=1 

<<7(110**11^ + |||V 2 d||Vd||| £ 2)||Vw||^2 (by(S2D) 

<\\\Q k u\\l* + ^Y^G{t) + C||Vu|| £4 . 

Next, we shall estimate the terms on the right-hand side of (14.171) . To bound the terms ||P|| £4 and ||Vn|| £ 4, 
it follows from Sobolev embedding, m , Holder’s inequality, (14.41) . (14.71) and (|4.18[) that 



ll^lll 4 + UVulli 4 <c(||VP|| 4 it + ||v 2 u|| 4 4 ) < C(\\eu\\\ + |||vd||v 2 d||| 4 4 ) 

<c||0||i a ||0>«||i a + q|vd||i a ||v 2 d||i4 

<C||^n|| 2 2 (F(t) + ||V«||l a ) + C||V 2 ci|| 2 2||V 3 d|| 2 2. (4.19) 


For the rest two terms ||V 2 d|| £4 and |||V<i| 2 || £ 4, by using Gagliardo-Nirenberg inequality, (14.41) and (14.181) . 
one has 


||V 2 d|| 4 £ 4 + |||Vd| 2 ||£ 4 <C||V 2 d|| 2 2||V 3 d|| 2 2 + C , ||Vd|| 2 2 ||V 2 d|| 2 2|||Vci||V 2 d||| 2 2 

<C||V 2 d|| 2 2(F(t) + ||Vu|| £4 ) + C'||V 2 d||| 2 ||V 3 d|| 2 2. (4.20) 

Hence, inserting (14.191) and (14.201) into (14.171) . one obtains 

^W + ^|V«|||2 + |||AVd||Vd||| 2 2 

<C{\\Q^u\\l* + ||V 2 d||i2)(P(t) + ||Vu||i2) + C||V 2 d||i2||V 3 d||i2. (4.21) 


Multiplying (14.211) by t 1 (i = 1,2), and then applying Gronwall’s inequality, it follows from (14.161) . (14.18[) . 
(14.21) . (14.41) and (14.41) that 


sup ifF(t)) + / f(||Vn || 2 2 + |||AVd||Vd||| 2 2 )dt 

0 <t<T Jo 


<C j t'-'F^dt + C [ iilV 2 d|| 2 21|V 3 d|| 2 2 dt + C f i\\e^u\\h + ^d^f^u^dt 
Jo Jo Jo 

<C [ T f-\\\giu\\l 2 + \\\V 2 d\\yd\\\l 2 )dt + C sup (t < - 1 ||Vw||| 2 ) / T ||V M || 2 2dt 
Jo 0 <t<T Jo 

+ C sup (p- 1 ||V 2 d|||2) / T t||V 3 d|| 2 2dt + C sup (tiVu|| 4 4 ) / T (||pia|| 2 2 + ||V 2 d|| 2 2)dt 
0 <t<T Jo 0 <t<T Jo 

<c f 2 +1|vd|| L 21|v 2 d|||21|v 3 d|| £2 )dt + c 

Jo 

<C f f- 1 (||pk|| 2 2 + ||V 3 d||i2)dt + C sup (t i_1 1| V 2 d|| 2 2 ) / ||V 2 d|| 2 2dt + C 

Jo 0 <t<T Jo 

<C. 


This together with (14.161) . (14. 18[) . (14.21) . (14.31) and (14.41) ensures (14.111) . 

Finally, it is easy to see that the estimate (14.111) combined with (14.71) implies (14.121) . This completes the 
proof of Lemma 14.31 □ 
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4.2 Higher order estimates 

In this subsection, we shall give some higher order norms estimates to the solution (g,u, P, d). To this end, 
we first list the following spatial weighted estimates on the density. 


Lemma 4.4 There exists a positive constant C depending on T such that 

sup \\gx a \\ L i < C(T). 

0 <t<T 


(4.22) 


Proof. Similarly to (13.1511 . we can prove (14.2211 with some suitable revisions, and we omit it. Moreover, we 
notice that (13.1411 still holds in R 2 . □ 


Lemma 4.5 There exists a positive constant C depending on T such that 

sup \\g\\mnw^+ l|V 2 U ||| 2 + ||VP || 2 2 + ||V 2 n||2f+||VP||2f + t(||V 2 n ||| 2ni9 + ||VP|| 2 2ni ,)) dt 


0 <t<T 


<C{T). 


(4.23) 


Proof. We first notice that for any r > 2, one follows from ED x and ED 4 that 

d 


dt 


livelier ^CIIVullLooUVellir. 


Using Gagliardo-Nirenberg inequality, and (l4~7l) . it follows that 


I|Vu||l- <C||Vn||^||V 2 U ||^ <C( ||HI 


2(4—1) 

L<* 


|||Vd||V 2 


iSr 17 


(4.24) 


(4.25) 


for all q £ (2, oo). On the one hand, by using Gagliardo-Nirenberg and Lemma T2.41 one has 


g(g-2) 


IIHU* <C\\gu \\ q L 2 1 \\gu\\ 

<c\\e^u\\ L i +c\\g?u\ 

From the above inequality, it is easy to see that 
rT / ,+i \ r T 


< c\\e*u\\£-'(\\e**\\L’ + ||VuM? 


g(g~ 2 ) 


2 g — 1 


g(g-2) 

UVullJ^, 


' 0 


\\gii\\ L q q +t\\gu\\ 2 L ^jdt <cj^ t ^ {t\\g^ u\\\^ ^ ^ (t\\^u\\h ) 2( * ^ dt + C J Q dt 

2q — l 


+ C J 0 + tWe^uWl 


<C sup 

0 <t<T 


(t\\giu\\l^j 


(29-1) -T 
29(9-1) / 


9 J + 9 Z -9-l 


f||Vii|| L 2+t 


dt 


+ C J (t||£> 2 u \\ 2 L 2 + t||Vu|| 2 2 ^ dt + C J (1 + \\g^u\\ 2 L2 )dt 


<C(T) 


(4.26) 


owing to (14.211 and (14.1111 . On the other hand, it follows from Holder’s and Gagliardo-Nirenberg inequalities, 
and (14.211 that 

' T / |||Vd||V 2 d|||)f + f|||Vd||V 2 d||| 2 , N ) dt 


-CJo +t(l|vd|| i 2 ,||v 2 d|| i2 ,) 2 ^dt 


<c 


l|Vd|II a ||V 2 d||L*||V 3 d|| ia ' 


Ml|Vd||« 3 ||V' J d|| La ||V d d|| La l dt 
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(4.27) 


<C ||V 3 d|| i2 ^ +t\\V 3 d\\ 2 L2 «)dt<C (||V 3 d|| 2 2 +f* + l)df <C(T). 


'0 \ /Jo 

Hence, combining (14.25[1 . (14.261) and (14.271) together, it follows that 

[ T ||Vu|| L =odt < C(T). 


(4.28) 


Thus, applying Gronwall’s inequality to (14.241) yields 


sup ||Vp|| L 2 nL9 < C(T). (4.29) 

0 <t<T 

Finally, it is easy to deduce from (14.71) . (14.261) . (14.271) . (14.21) and (14.41) that 

r T / 9 +i 9 ±i \ 

J o (j|V 2 u||£ 2 + ||VP||£ 2 + ||V 2 u||J +||VP|| L J +t(\\V\\\l 2nLq + \\VP\\l 2nLq )j dt < C(T). 

This together with (13.61) and (14.291) ensures (14.231) . This completes the proof of Lemma l4~5l □ 

In the following Lemma, we shall give some spatial estimate on Vp, Vd and V 2 d, which are crucial to 
derive the estimates on the gradients of both ut , d* and Vd(. 


Lemma 4.6 There exists a positive constant C depending on T such that 


sup \\gx a \\ L ^nH^nw l ’i < C(T), 

0 <t<T 


sup ||Vdx^ || 2 2 + [ \\V 2 dx%\\ 2 L 2 dt<C(T), 

0 <t<T Jo 


(4.30) 

(4.31) 


and 


sup t|| V 2 da; 2 1|| 2 + f f||V 3 da; 2 || 2 2 dt < C(T). (4.32) 

o <t<T Jo 


Proof. With (14.221) in hand, the proof of (14.301) is exactly the same as [351 Lemma 3.6]. Similarly to (13.51) . 
we can prove (14.311) with some suitable revisions, and we omit it for simplicity. 

Now, multiplying (13.81) by A Vdx a and integrating by parts, it follows that 


1| V 2 dS^ || 2 2 +|| V 3 dx^ || 2 2 = -J\/d t W 2 dVx a dx + jV{u ■ Vd)VAdx a dx 

- j V(|Vd| 2 d)VAdi a da; -2 J V 3 dV 2 dVx a dx - J \V 2 d\ 2 V 2 x a dx 
<i|| V 3 dS^ ||| 2 + J |AVd||V 2 d|WMx+ J |u||Vd||V 3 d|V:r a dx+ J \u\\S7d\\V 2 d\V 2 x a dx 
+ J\u\\\7 2 d\ 2 Vx a dx + J |Vd| 3 |V 2 d|Va; Q d:r + J \Vd\\V 2 d\ 2 Vx a dx 
+ J \Vu\ 2 \Vd\ 2 x a dx + J |Vd| 6 5 0 dx + J \Vd\ 2 \V 2 d\ 2 x a dx + J \V 2 d\ 2 V 2 x a dx 


aI|| V 3 dit ||| 2 + Kl + K 2 + ---+ K w . 

We estimate all the terms I\i(i = 1, 2, • • ■ , 10) on the right side of (14.331) term by term as follows, 


(4.33) 


A'i < J\VAd\\V 2 d\x a dx < ^|| V 3 di^ ||| 2 + C\\S7 2 dx^ ||| 2 , 

K 2 < J |u||Vd||V 3 d|x a -tdx< ^||V 3 di^||i 2 + C7|| V 2 dit||2 4 || ui -|||2 a 

<^!|V 3 dxt|| 2 2 +C7||V 2 d5;t|| i2 (||V 3 dxf|| L2 + ||V 2 dVst|| i2 )(||^ u ||2 2 + ||V M || 2 2 ) 
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<i||V 3 ^t||2 2+C || V 2^f||2 2) 

IC 3 <C\\V 2 dx^\\ L 2\\S7d,x^\\ L i\\ux~^\\ L i 

<||v 2 d®^|U a ||Vd®’|| La (||v 2 d®^|| ia + HVdVx^xil^lU 2 + iivuiuo 

<C(l + ||V 2 d^||| 2 ) 

K 4 <C\\V 2 dx* || L 2 ||V 2 d.T 2 || L 4||'uS _ 3||l4 

<||VW||| 2 ||V 2 ^t||J 2 (||V 3 dxf || L 2 + ||V 2 dVxt|| i2 )i(||^ u || L2 + IIVulUO 
<^||V 3 ^f|| 2 2 +C||V 2 dxt|| 2 2) 

K S <C J |Vd| 3 |V 2 d|i°da; < C||V 2 di^|| L2 ||Vc?x 2 IM|Vd||| 8 

<C\\S7 2 dx% \\ L 2 || V<£e^ ||| 2 (|| V 2 dx^ \\ L 2 + || VdVa;^ ||l 2 )^ II Vd||| 2 || V 2 d||f 2 
<C(l + ||V 2 ^f|| 2 2 ), 

K e <cj |Vd||V 2 ci|Vx a da;<q|VdS^|| L 2||V 2 dx®|||4<C'||V 2 dx^|| L 2(||V 3 dxt|| L2 + ||v 2 d5^|| L 2) 

<l||V 3 ^f||2 2+C || V 2 ^f||2 2) 

K r <C\\Vdx* ||i 4 ||Vu||| 4 < C||Vctet|| L2 (||V 2 (&5|| L 2 + ||VdVS5|| L 2)||Vw|| L2 ||V 2 u|| i 2 
<C|| V 2 m || 2 2 + C(1 + \\S7 2 dx^ \\l 2 ), 

Ks + I< 9 <C||Vd||lsIIV 2 c&t ||2 4 + C||Vd||£ 4 ||V 2 d®S ||2 4 

<C(||Vd|| L 2 ||V 2 d||| 2 + ||Vd|| L 2 ||V 2 d|| L2 )||V 2 dxt|| i 2 (||V 3 dxt|| L2 + ||v 2 dVst|| i 2 ) 
<l||V 3 dxf|| 2 2 +C||V 2 dxf|| 2 2) 
ivio < J \V 2 d\ 2 x a dx <C\\V 2 dx%\\ 2 L2 , 

where we have used Holder’s and Gagliardo-Nirenberg inequalities, (14.21) . (14.41) . (13.141) and (14.311) . Inserting 
these estimates Ki(i = 1,2, • ■ • , 10) above into (14.331) . after by using (14.211) . we have 

^l|V 2 d* f ||| 2 -H|V 3 d55||| 2 < C\\V 2 dx% ||| 2 + C(1 + ||V 2 U ||i 2 ). 

Multiplying the above inequality with t, and then using Gronwall’s inequality and (14.23|) . we prove (14.321) . 
this completes the proof of Lemma 14.61 □ 

Lemma 4.7 There exists a positive constant C depending on T such that 

sup t (\\^u t \\ 2 L2 + \\dt \\ 2 H1 + IIv 3 d||| 2 )+ [ T t(\\Vu t \\ 2 L 2 + HVdtll^) dt < C(T). (4.34) 

0 <t<T v ' J 0 

Proof. We shall first prove that 

£ (We'utWb + WVdtWh) d t < C(T). (4.35) 

Indeed, it is easy to see 

W^ut || 2 a <|[^«|| 2 3 + ||e*M|Vu||| 2 2 < \\g^u\\ 2 L2 +C\\g^u\\ 2 L6 \\Vu\\ 2 L3 
<\\qH\\ 2 l2 +C(|| e * u || 2 a + ||V M ||i 2 )||V u ||| 2 ||V 2 w ||| 2 

<|[e**ll£ 2 +C'(l + ||V 2 u|| 2 3 ) (4.36) 
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due to Lemma [2~4l m and (14.41) . By using (13.81) and Lemma [2751 it follows that 


||Vd t ||£ 2 <C(||V 3 d || 2 


L 2 


<C(\\V 3 d\\i 2 + 

73 J||2 


<C'(||V 3 d || i 2 
<C(\\V 3 d\\ 2 L2 
<C(\\V 3 d\\ 2 L2 


l|V«||Vd||| 


L 2 


H|v 2 d|||i 2 + |||Vd| 


3| li 2 


|||Vd||V 2 d||| 


L 2 ' 


3 ^'" 2 1 lVu||| 4 ||Vd ||| 4 + \\ux-l\\ 2 L s\\V 2 dx%\\ L2 \\V 2 d\\ Li + ||Vd||£a + |||Vd||V 2 d||| 2 2 ) 

I V 2 m || 2 2 + ||V 2 ^t||| 2 + ||nx-t|| 4 i8 ||V 2 d|| 2 4 + 1) 


V 2 «|| 2 2 + II V 2 dx%W 2 L2 + (||^u ||| 2 + II Vm|| 2 2 ) 2 || V 2 d||^4 + 1) 
V 2 m || 2 2 + 1) 


(4.37) 


owing to (14.21) . (14.41) and (14.311) . Combining (14.361) . (14.371) . (14.21) and (14.23|) together, it is easy to see that 
(14.351) holds. 

Now, multiplying (13.251) by u t and integrating the resulting equality by parts over R , we obtain after 
using tEDl and the divergence free condition ( 11 . 11 ) t that 


2 dt' 


Q*ut\\l 2 + l|Vut||| a < Cj g\u\\u t \ (|Vw t | + |Vu | 2 + |u||V 2 M |) dx 

+ C j £>|n| 2 |Vu||Vu t |d:r + C J p|wi| 2 |Vu|dx + C J |Vd||Vd t ||Vu t |da: 


(4.38) 


=Mi + M 2 + M 3 + M 4 . 

By using (14.21) . (14.41) . Holder’s and Gagliardo-Nirenberg inequalities, it follows that 

Ml <C'|| e 5 W || L6 || e ,5 Mt || 2 2 ||^ Ut || 2 6 (||Vu t || i 2 + ||Vu||| 4 ) +C||p2M|||i 2 ||£»5 Ut || 2 2 || £) 5 Ut || 2 6 ||V 2 u|| i 2 

<C’(l + ||Vu|| 2 2 )||pii tt ||| 2 (||^ Ut || i2 + ||v ltt || i 2 )^(||V Ut || i 2 + ||V U || 2 2 + ||Vn|| i 2 ||V 2 U || i 2 + ||V 2 n|| i 2 ) 

<g||Vut ||| 2 + C( 1 + ||£> 2 Wi ||| 2 + ||V 2 u|| 2 2 ), 

M 2 +M 3 < C’||^ u || 2 a||V M || i4 ||Vn t || i2 + C\\^u t \\l 6 \\g^u t \\l 2 \\Vu\\ L , 

<g||VM t || 2 2 + C( 1 + ||£> 2 Ut ||| 2 + II V 2 u||| 2 ), 

M 4 <C||Vd|| L 4||Vd t |U 4 ||Vu t || L 2 < i||V«t||l 2 +C||Vd|| i 2||V 2 d|| L 2||Vdt|| L2 ||V 2 d t || L 2 

u 

where the positive constant C 3 is defined in the following (14.401) and (14.431) . Substituting the estimates of 
Mi{i = 1, 2, • • • ,4) into (14.381) . it follows that 


^Il0 2 u t|]l 2 + l|Vwt ||| 2 < C'(||p 2 M t ||| 2 + ||Vdt||| 2 ) + + ^ 


l|V 2 d t |li 2 +C(||V 2 n||i 2 + l). 


(4.39) 


Differentiating (ED3 with respect to time variable t, multiplying the resulting equality with dt and then 
integrating by parts over R , we have 

\^\\dt\\h + \\Vdt\\l* < c J l u t||Vd||df |dx + C j |Vdt||Vd||dt|da; + C J |Vd| 2 |d t | 2 dx 

—M§ + AIq + M 7 . 


By using Holder’s and Gagliardo-Nirenberg inequalities, (14.21) . (14.41) . (13.141) and (14.311) . we have 

Me <\\utx~^\\ L i\\S7dx^\\ L 2\\dt\\ L i < ||l 2 + ||VM t || L 2 )||Vd5^|| L2 ||d t || 2 2 ||Vdt||2 2 

< i l|V Wt ||| 2 + C{\\q 2 u t \\\ 2 + ||dt ||| 2 + ||Vd t ||| 2 ), 

Me+M 7 <^\\Vd t \\ L , + C\\Vd\\U\d t \\ 2 Li < i||Vd t || i 2 +G||Vd|| i 2 ||V 2 d|| i 2 ||d t || i 2 ||Vdi || i2 
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Hence 


<^d t \\ L 2+C\\d t \ 


2 

L 2 * 


^IMtlli 2 + HVdJl, <C 3 \\Vu t \\ 2 L 2 + C(||£. = u t ||| 2 + ||d t |||2 + ||Vdt|||2). (4.40) 

Differentiating (13.81) with respect to time variable t ensures 

Vd tt - AVd t = —V(u • Vd)t + V(|Vd| 2 d) t . (4.41) 

Multiplying (14.411) by Vd t , and integrating the resulting equality over K 2 , it follows that 

— WVdtWh + W^dtWbKC j |Vu t ||Vd||Vd t |dx + C J \Vu\\Vd t \ 2 dx 

+ C J \u t \\Vd\\V 2 d t \dx + C J \Vd\ 2 \d t \\V 2 d t \dx + C J |Vd||Vd t ||V 2 dt|da; 

— Ms + Mg + Mio + Mu + M 12 . (4.42) 


By using Holder’s and Gagliardo-Nirenberg inequalities, (14.21) . (14.41) . (13.141) and (14.311) we have 
M 8 <\\Vu t \\ L *\\Vd t \\ L 4Vd\\ L * < ^||Vu t || 2 2 +C||Vd|| L 2 ||V 2 d|| i 2 ||Vd t || i 2 ||V 2 d t || L 2 

<|l|V^||i 2 + |||V 2 d t ||| 2 + C'HV^Hi,, 

M 9 <C\\Vu\\ L 2 \\Vd t \\ 2 Li < C'||Vu|| L 2 ||Vd t || L 2 ||V 2 d t || i 2 < i||V 2 d t || L 2 + C\\Vd t \\ L 2 , 

Mio <C'y’|u t x- 22 f i ||Vdx^| 2 ^ i |Vd|^|V 2 d t |da ; <C , ||« t x- 2s f i || i8 .||Vdxt|| J i ||Vd||^||V 2 d||^||V 2 d t || i 2 
<gll^ 2 ^ll |2 + C\\u t X 4 |||sa < -|!V 2 d t ||| 2 + C(\\QlU t \\ L 2 + ||Vu t ||| 2 ), 

Mn <i||V 2 d*|| 2 2 +C||Vd|| 2 8 ||d t || 2 4 < i||V 2 d t || 2 2 +C'||Vd||i 2 ||V 2 d||| 2 ||d t || i 2||Vd t || i 2 


<^' 2 dt\\l,+C(\\d t \\ 2 L2 + \\Wd t \\ 2 L2 ), 

-Mi 2 <^l|V 2 d t || 2 2 +C||Vd|| 2 4 ||Vd t || 2 4 < i||V 2 d t || 2 2 +C'||Vd t || 2 2. 

Inserting the estimates of M t (i = 8,9, • ■ ■ , 12) into (14.421) . it follows that 


^||Vd t || 2 2 + ||V“d t || 2 2 < C , 3(||e 2 ut|| i 2 + ||Vu t || 2 2 ) + C(||d t ||| 2 + ll^tlll 2 )- (4.43) 

Next, multiplying (14.391) by 2 C 3 + 1 and adding the resulting inequality with (14.401) and (14.431) . we have 

— ((2C 2 + l)||^w t |li 2 + \\d t \\ 2 H i) + ||V Ut || La + -||Vd t ||ffi 
<C{ 1 + WgiutWh + l|Vd t || 2 ffl ) + C{1 + \\V 2 u\\ 2 l2 ). 


Multiplying the above inequality by t, we obtain 

,.T 

sup t(\\ Q 1 2u t \\ 2 L2 + \\X7d t \\ 2 H i)+ / t (\\Vu t \\ 2 L2 + \\Wd t \\ 2 Hl ) dt < C{T), (4.44) 

0 <t<T v 7 JO 

after by using Gronwall’s inequality and (14.81) . 

Finally, it follows from Holder’s and Gagliardo-Nirenberg inequalities, (14.21) . (14.41) and (13.141) that 

||V 3 d|| 2 2 <C(||Vd t |li2 + |||Vu||Vd||| 2 2 + HM|V 2 d||| 2 2 + |||Vd| 3 || 2 2 + |||V 2 d||Vd||| 2 2) 
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<C(\\S7d t f L2 + ||Vn|| 2 4 ||Vd|| 2 4 + ||t t ®-S||i.||V 2 d®S|| ia ||v 2 d|| L 4 + ||Vd||£. + |||Vd||V 2 d|||| 2 ) 
<C(||Vd t ||| 2 + ||V 2 M ||| 2 + \\V 2 dx%\\li + \\ux~? Hls||V 2 d||i 4 + ||Vd||i 2 ||V 2 d|| 4 i2 + ||V 2 d|||3||Vd||ia) 
<G(||Vd t || 2 2 + ||V 2 u|l 2 2 + ||V 2 ^^|| 2 L2 + (||^ U || 2 L2 + ||Vn||| 2 ) 2 ||V 2 d|| L2 ||V 3 d|U 2 
+ ||V 2 d||l 2 + ||V 2 d|| i2 ||V 3 d|| i2 ) 

<\\\^d\\\ 2 + C(||Vd t || 2 2 + ||V 2 u|| 2 2 + ||V 2 d|| 4 2 + ||V 2 da:t ||| 2 + ||V 2 d|| 2 2 ) 

<^l|v 3 d||! 2 + C(||vd t || 2 2 + ||V 2 U || 2 2 + ||v 2 ctef I|| 2 + ||v 2 d|| 2 2 ), 

which combined with (14.31) . (14.121) . (14.321) and (14.441) indicates (14.341) . This completes the proof of Lennna l4.7l 

□ 

4.3 The proof of Theorem 11.21 

In this subsection, with Theorem 13.11 and the a priori estimates obtained in Subsections 4.1 and 4.2 at hand, 
we shall give the proof of Theorem 11.21 

By Theorem l3.ll we know there exists a T* > 0 such that the Cauchy problem of system (11.11) (11.21) admits 
a unique strong solution (g, u, P, d) on R 2 x (0, T*]. In what follows, we shall extend the local solution to all 
the time. 

Set 


T* = sup { T\(g , u, P , d) is a strong solution to (11.11) (11.21) on R 2 x (0, T]} . (4.45) 

First, for any 0 < t <T* <T <T* with T finite, one deduces from (14.21) . (14.41) . (14.121) and (14.341) that for all 

q >2, 

Vu, Vd, V 2 d G C ([r, T]; L 2 (R 2 ) n L^R 2 )), (4.46) 

where one has used the standard embedding 

L 00 (T,T;P 1 (R 2 ))nP 1 (r,r;P- 1 (R 2 )) C(r, T; L 9 (R 2 )) for all q G [2,oo). 

Moreover, it follows from (14.231) . (14.301) and [32) Lemma 2.3] that 

e g c([o, t] ; l 3 (r 2 ) n h\ r 2 ) n w li9 (R 2 )). (4.47) 

Now, we claim that 


T* = 00. 

Otherwise, if T* < oo, it follows from (14.461) . (14.471) . (14.21) . (14.41) . (14.301) and (14.311) that 

{g,u, P, d)(x,T*) = lim (g,u, P,d)(x,t) 

t^T* 

satisfies EzD at t = T*. Moreover, using m and (13.61) with p = 1, it follows that 

f g(x, T*)dx = f go{x)dx = 1. 

Jr 2 Jr 2 

Notice that there exists N 0 > 0, it is easy to see that 

[ g{x,T*)dx > \ f g(x,T*)dx > 

Jr 2 1 Jb Nq z 


(4.48) 


33 




























Thus, we can take (g, u , P, d)(x , T*) as the initial data, and Theorem l3 .1 l implies that one could extend the local 
solutions beyond T*. This contradicts the assumption of T* in (14.4511 . Hence, we prove (14.481) . Furthermore, 
from (14.31) . (14.111) . (14.121) and (14.341) . one obtains that (11.101) holds. This completes the proof of Theorem ll.2l □ 
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